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ABSTRACT 

We present a catalog of 1439 young stellar objects (YSOs) spanning the 1.42 deg 2 field surveyed 
by the Chandra Carina Complex Project (CCCP), which includes the major ionizing clusters and the 
most active sites of ongoing star formation within the Great Nebula in Carina. Candidate YSOs were 
identified via infrared (IR) excess emission from dusty circumstellar disks and envelopes, using data 
from the Spitzer Space Telescope (the Vela-Carina survey) and the Two-Micron All Sky Survey. We 
model the 1-24 /im IR spectral energy distributions of the YSOs to constrain physical properties. Our 
Pan-Carina YSO Catalog (PCYC) is dominated by intermediate-mass (2 M Q < m < 10 M Q ) objects 
with disks, including Hcrbig Ac/Be stars and their less evolved progenitors. The PCYC provides a 
valuable complementary dataset to the CCCP X-ray source catalogs, identifying 1029 YSOs in Carina 
with no X-ray detection. We also catalog 410 YSOs with X-ray counterparts, including 62 candidate 
protostars. Candidate protostars with X-ray detections tend to be more evolved than those without. 
In most cases, X-ray emission apparently originating from intermediate-mass, disk-dominated YSOs is 
consistent with the presence of low-mass companions, but we also find that X-ray emission correlates 
with cooler stellar photospheres and higher disk masses. We suggest that intermediate-mass YSOs 
produce X-rays during their early pre-main sequence evolution, perhaps driven by magnetic dynamo 
activity during the convective atmosphere phase, but this emission dies off as the stars approach the 
main sequence. Extrapolating over the stellar initial mass function scaled to the PCYC population, 
we predict a total population of >2 x 10 4 YSOs and a present-day star formation rate (SFR) of 
>0.008 M Q yr _1 . The global SFR in the Carina Nebula, averaged over the past ~5 Myr, has been 
approximately constant. 

Subject headings: circumstellar matter — ISM: individual (NGC 3372) — stars: formation — stars: 
luminosity function, mass function — stars: pre-main sequence — X-rays: stars 



1. INTRODUCTION 

The Great Nebula in Carina (NGC 3372) is among 
the most active massive star-forming complexes in the 
Galaxy. Over 70 O-type stars ionize a giant H II region 
(Smith 2006a) and power a bipolar superbubble that 
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has erupted from a 100-pc scale giant molecular cloud 
complex (Smith ct al. 2000). The Carina Nebula is often 
considered to be the nearest analog of extragalactic star- 
burst regions (d — 2.3 kpc; Smith 2006b), and it is also 
the least obscured. Carina thus presents a rich target for 
the study of resolved young stellar populations evolving 
in an environment dominated by some of the most mas- 
sive stars known. While numerous studies have targeted 
individual rich ionizing clusters in Carina (Feinstein et al. 
1980; Massey & Johnson 1993; DeGioia-Eastwood et al. 
2001; Ascenso et al. 2007; Sana et al. 2010), global obser- 
vations of the stellar populations across the complex are 
challenging. Due to its large angular size (>2° x 1.5°), 
it is difficult to achieve the requisite wide observational 
coverage of the Carina Nebula at high spatial resolu- 
tion. Its position near the Galactic midplane (6 = — 0?6) 
means that wide-field images at visual and infrared (IR) 
wavelengths are dominated by contaminating popula- 
tions of unassociated stars. The Chandra Carina Com- 
plex Project (CCCP), a 1.2 Ms total X-ray integration 
covering a 1.42 deg 2 field in the center of the Carina 
Nebula, provides an unprecedented opportunity to over- 
come these observational challenges and study the young 
stellar populations from a global perspective, yet at the 
detailed level of resolving individual stars (Townsley et 
al. 2011). 

X-ray observations of star-forming regions are highly 
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sensitive to low-mass pre-main-sequence (PMS) stars 
(T Tauri stars), which generate hard X-rays via 
convectively-driven magnetic flaring activity (Preibisch 
et al. 2005) and to OB stars that produce soft X-rays pri- 
marily through micro-shocks in their strong stellar winds 
(Lucy & White 1980; Feldmeier et al. 1997). Both clas- 
sical T Tauri stars with circumstellar disks and proto- 
stars with accreting envelopes are observed to be sig- 
nificantly less luminous in X-rays compared to diskless, 
weak-lined T Tauri stars (Stassun et al. 2004; Telleschi 
et al. 2007; Prisinzano et al. 2008). The CCCP observa- 
tions are expected to be insensitive to intermediate-mass 
(2 M Q < m < 10 M ) young stellar objects (YSOs), 
including Hcrbig Ac/Be stars with disks and their less 
evolved progenitors, which lack a known mechanism for 
the production of strong X-ray emission (Stelzer et al. 
2006). This important stellar population traces recent 
and ongoing star formation in massive complexes like 
Carina. Because such stars have high bolomctric lu- 
minosities and possess dusty disks and envelopes that 
reprocess the stellar radiation into strong mid-IR emis- 
sion, surveys with the Spitzer Space Telescope are very 
sensitive to them (Koenig et al. 2008; Robitaille et al. 
2008; Kang et al. 2009; Povich & Whitney 2010, here- 
after PW10). The entire Carina Nebula was included in 
the Spitzer Vela-Carina survey. This complementarity 
between the IR and X-ray observations motivates us to 
search for candidate YSOs throughout the Carina Neb- 
ula using Vela-Carina survey data in conjunction with 
supplementary IR datasets. 

The Carina complex contains at least 6.7 x 10 5 M 
in molecular gas mass (Grabelsky et al. 1988) but hosts 
few bright embedded IR sources or other signposts of 
ongoing massive star formation; hence historically the 
Carina Nebula has been regarded as a relatively evolved 
H II region, with a current star formation rate (SFR) 
that is negligible in comparison with its historical rate 
(Smith & Brooks 2008, and references therein). Within 
the past decade, a new view of star formation in Ca- 
rina has emerged, in which active star formation has mi- 
grated to the outer reaches of the nebula, in particular 
the region known as the South Pillars, where remnant 
molecular cloud material has been shredded and sculpted 
into myriad clumps, pillars, and elephant-trunk struc- 
tures by stellar winds and radiation, the two dominant 
feedback mechanisms produced by massive stars during 
their lifetimes (Megeath et al. 1996; Smith et al. 2000; 
Rathborne et al. 2004; Smith & Brooks 2007). Smith 
et al. (2010b, hereafter SP10) presented the first Spitzer 
imaging study of Carina using the Infrared Array Cam- 
era (IRAC; Fazio et al. 2004) and found over 900 candi- 
date YSOs widely distributed throughout two large fields 
covering the South Pillars and a portion of the opposite 
wall of the southern supcrbubblc lobe. Carina remains a 
highly active star-forming region, but it may be unusual 
among Galactic giant molecular glouds (GMCs) in the 
sense that the dominant mode of current star formation 
activity appears to be triggering by feedback from the 
local massive stellar population (SP10). 

Our present study expands upon the Spitzer point- 
source results presented by SP10 in several ways: (1) The 
Vela-Carina observations extend Spitzer coverage to the 
entire Carina Nebula, including the central ionizing clus- 
ters and regions to the North. The wide-field coverage 



enables us to analyze relatively empty "control" fields 
far from the Carina GMC and correct for contamination 
in our YSO sample from unassociated sources with sim- 
ilar mid-IR colors. (2) We incorporate 24 /zm photome- 
try from the Multiband Imaging Photometer for Spitzer 
(MIPS; Ricke et al. 2004) to constrain the evolutionary 
stages and luminosities of individual YSOs better than 
what is possible with IRAC photometry alone. (3) We 
analyze the subsample of YSOs that match CCCP cata- 
log sources to investigate the physical properties of can- 
didate X-ray-emitting protostars and intermediate-mass 
PMS stars with disks. 

In this paper, we present a highly-reliable catalog of 
1,439 predominantly intermediate-mass YSOs through- 
out the Carina Nebula, selected via their IR excess emis- 
sion. This Pan-Carina YSO Catalog (PCYC) comple- 
ments and augments the CCCP X-ray catalog (Broos et 
al. 2011a). The data analysis steps used to construct 
the PCYC are summarized in §2, and the methodology 
used to model the physical properties of the YSOs is de- 
scribed in §3. In §4, we present the ensemble properties 
of the YSO population, in §5 we analyze the subset of 
our sample associated with X-ray emission, and in §6 we 
highlight five noteworthy YSO clusters and sub-groups 
and several interesting individual YSOs. We discuss our 
results in the context of the star formation history of the 
Carina Nebula and the physical mechanism(s) powering 
X-ray emission from YSOs in §7. We summarize our 
conclusions in §8. 

2. DATA ANALYSIS AND YSO SAMPLE SELECTION 

2.1. Observations 

The Vela-Carina survey (Spitzer Proposal ID 40791, PI 
S. R. Majewski), was modeled after the Galactic Legacy 
Infrared Mid-Plane Survey Extraordinaire (GLIMPSE; 
Benjamin et al. 2003). It used the 4 mid-IR bands of 
Spitzer/IRAC (3.6, 4.5, 5.8, and 8.0 /xm) to image 86 deg 2 
of the Galactic plane in a strip 2° high in latitude over 
longitudes 255° < I < 295° (Zasowski et al. 2009), an 
area than included the Carina complex. PSF-fitting pho- 
tometry was performed using the GLIMPSE point-source 
extractor, a modified version of DAOPHOT (Stetson 
1987) optimized for crowded fields with strong, varying 
background nebular emission. The GLIMPSE extractor 
produces two source lists, the 99.5%-reliable Point Source 
Catalog and the more complete Point Source Archive. 14 
The depth of the Vela-Carina photometry is very simi- 
lar to that of GLIMPSE, reaching [3.6] <15.5 mag (less 
sensitive in regions of bright nebular emission). Both 
the Catalog and the Archive sources have been matched 
(bandmcrged) with JHKs photometry from the Two- 
Micron All-Sky Survey Point Source Catalog (2MASS; 
Skrutskic et al. 2006), and the resulting source lists pro- 
vide photometry in 7 IR bands, from 1.2 /im through 
8.0 /im. We use primarily the highly-reliable Catalog for 
this work, as experience shows that this drastically re- 
duces the number of false-positives, candidate YSOs that 
are generated as the result of systematic photometric 
errors. The Vela-Carina Catalog contains 60,515 point 
sources inside the CCCP survey area (Table 1). We do re- 
trieve a few sources from the Archive that are matched to 

14 For details of the data processing and products, go to 
http://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/doc/glimpsel_dataproc 
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bright candidate YSOs detected by the Midcourse Space 
Experiment {MSX), as described in §2.4 below. 

The CCCP catalog contains over 14,000 point sources 
detected by the Advanced CCD Imaging Spectrometer 
(ACIS; Garmire et al. 2003) of the Chandra X-ray Ob- 
servatory. Broos et al. (2011a) matched the Vela-Carina 
Archive to the CCCP X-ray catalog and found 6,543 
counterparts to CCCP sources. The Vela-Carina Cata- 
log is a subset of the Archive, and provides mid-IR coun- 
terparts for 4,664 CCCP sources (Table 1). 

The Carina Nebula was observed with Spitzer/MIPS 
during late July 2007 as part of program GO-30848 
(MIPSCAR; P.I. N. Smith). The program consisted of a 
series of MIPS scans that covered most of the area of the 
star forming region but avoided the bright emission of 
the Homunculus nebula around r] Carinae. This yielded 
a large mosaic image of the 24 /j,m emission that contains 
the entire CCCP survey area. 

2.2. Identification of Candidate YSOs 

YSOs possess dusty circumstellar disks and infalling 
envelopes that reprocess radiation from the central stars, 
producing characteristic IR excess emission. YSOs can 
be identified in broad-band photometric imaging obser- 
vations via their IR colors or spectral energy distribu- 
tions (SEDs). Our procedure for identifying candidate 
YSOs is adapted from Povich et al. (2009), SP10, and 
PW10. Here we outline the main steps used in con- 
structing our highly-reliable catalog of candidate YSOs, 
and we summarize the number of sources found within 
the CCCP survey area for each step in Table 1. Us- 
ing a x 2 -minimization SED fitting tool (Robitaille et al. 
2007), we fit reddened Castelli & Kurucz (2004) stel- 
lar atmospheres (allowing Ay to vary between and 
40 mag according to the extinction law of Indcbetouw 
et al. 2005) to the SEDs of 210,742 sources in the Vela- 
Carina Catalog located within a Al x Ab = 2? 8 x 2° field, 
54,155 of which fell within the CCCP survey area (Ta- 
ble l). 15 Only sources detected in Adata > 4 of the 7 
Catalog bands (A = 1-8 /im) were fit. To avoid biasing 
the fits for sources where the photometric uncertainties 
may have been underestimated in the Catalog, before fit- 
ting any models we conservatively reset the uncertainties 
to a minimum value of 10%. We considered sources for 
which Xo/-^data < 2 (x§ is the unreduced goodness-of-fit 
parameter of the best-fit model; SP10) to be well-fit by 
stellar atmosphere SEDs; among sources in the CCCP 
field that were fit with SED models, 93.4% were con- 
sistent with normally-reddened stellar photospheres (Ta- 
ble 1) and hence were removed from consideration for the 
PCYC. 

We then filtered out sources with "marginal" IR ex- 
cess emission, where the excess appears in only the single 
IRAC [5.8] or [8.0] band, using the SP10 color criteria, 
modified to de-redden the [3.6] — [4.5] colors of back- 
ground stars viewed through the Carina cloud by up to 
Ay = 20 mag. Using [A] to denote magnitudes in the var- 
ious IRAC bands and 6([\i] — [A.,-]) for the uncertainties 
on the colors computed from the (minimum 10%) uncer- 
tainties on Catalog flux densities, the selection criteria 

15 The full field analyzed is much larger than the CCCP field; 
regions of active star formation identified in this wider field will be 
the subject of a future paper. 



were as follows: 

[3.6] - [4.5] > £([3.6] - [4.5]) + £([3.6] - [4.5]) 
OR | [4.5] - [5.8]|><5([4.5] - [5.8]) 
AND [5.8] - [8.0] ><J([5.8] - [8.0]). 

The color excess used for the de-reddening was calculated 

as 

£([3.6] - [4.5]) - A v ( " 3 - 6 ~ K4 - 5) = 0.0135A y , 

Ky 

where the k\ are opacities given by the extinction law 
(Indcbetouw et al. 2005). This step is crucial to dis- 
tinguish intrinsically red objects from sources affected 
by potential systematic photometric errors (Povich et al. 
2009); IRAC [5.8] and [8.0] are less sensitive to point 
sources and more affected by nebular emission compared 
to IRAC [3.6] and [4.5]. We discarded 1,923 of these 
"marginal" IR excess sources, leaving 1,646 "reliable" 
intrinsic IR excess sources, which form the basis for our 
sample of candidate YSOs (Table 1). Throwing out the 
majority of the possible IR excess sources in this manner 
may seem overly conservative, but we note that among 
the CCCP-matched subsamplc, only 1/3 of sources were 
caught by this filter (Table 1). The Vela-Carina survey 
images are dominated by K- and M-type giants in the 
Galactic field that are not expected to emit X-rays (Ros- 
ner et al. 1995). If the majority of marginal IR excess 
sources were actually due to emission from disks, perhaps 
transitional disks with large inner holes, we would expect 
the fraction of sources rejected at this step among the 
X-ray detected subsample to be similar to, if not higher 
than, the fraction rejected from the full sample. Instead, 
the significantly lower frequency of marginal IR excess 
sources with X-ray counterparts affirms that the major- 
ity of these sources are unassociated field stars. The 213 
CCCP-matched sources in Table 1 rejected on account 
of marginal IR excess are examined again by Povich et 
al. (2011, hereafter Pll). 

2.3. MIPS 24 Aperture Photometry 

Our data analysis is tailored to the goal of fitting 
YSO model SEDs from Robitaille et al. (2006, hereafter 
RW06) to the available IR photometry of each candidate 
YSO. Because the SEDs of highly embedded YSOs peak 
in the thermal IR near 100 fim, photometry at A > 10 /im 
is often required to constrain bolometric luminosities and 
to distinguish envelope-dominated (Stage 0/1) from disk- 
dominated (Stage II) YSOs (RW06; Indebetouw et al. 
2007). Following PW10, we measured MIPS 24 ^m flux 
densities using aperture photometry. We located the po- 
sition of each candidate YSO and centered an extraction 
aperture of radius 3.5" and background annulus of inner 
and outer radii 7" and 13", respectively, on the MIPS 
24 /zm mosaic of Carina. We extracted 24 /xm flux den- 
sities for each source, estimating the background level 
using the Daophot MMM algorithm (Stetson 1987). Our 
choice of extraction aperture required an aperture correc- 
tion 16 of 2.8. The uncertainty introduced by the aperture 
correction is more than offset by the greater accuracy of 

16 MIPS Instrument Handbook vl.0, 

http://ssc.spitzer.caltech.edu/mips/mipsinstrumenthandbook/ 
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Table 1 

Summary of Vela-Carina Point Sources in CCCP Field 



Sources 



All 



CCCP-matchcd 



Description 



In Vela-Carina Catalog 60,515 

Fit with SED Models 54,155 

Stellar Photospheres 50,586 

Possible IR Excess 3,569 

Marginal IR Excess 1,923 

Reliable IR Excess 1,646 

In PCYC a M39" 

Stage 0/1 247 

Stage II 478 

Stage III 83 

Ambiguous 631 



4,664 Detected IRAC point sources, 99.5% reliable 

4,080 Detected in >4 of 7 2MASS+IRAC bands 

3,444 Well-fit by stellar atmosphere SEDs 

636 Poorly-fit by stellar atmosphere SEDs 

213 Excess in only IRAC [5.8] or [8.0] band 

423 Intrinsically red objects 

410 YSO catalog, corrected lor contamination 

62 Dominant SED models: infalling envelopes 

190 Dominant SED models: optically thick disks 

13 Dominant SED models: optically thin disks 

144 Inconclusive SED model classification 



Includes 3 Vela-Carina Archive sources with MSX counterparts. 



the local background determination and the ability to 
separate close sources (PW10). 

PW10 analyzed MIPS GAL 24 /im survey data (Carey 
ct al. 2009) covering the M17 SWex IR dark cloud. The 
MIPSCAR observations are shallower than the MIPS- 
GAL observations and affected by extremely bright neb- 
ular emission produced by warm dust in the Carina 
H II region. While ~65% of the YSOs cataloged by 
PW10 were detected at 24 /jm, the detection fraction 
is only ~10% for the PCYC. We can, however, derive 
useful upper limits on 24 /zm flux densities in the ma- 
jority of cases. Upper limits provide critical constraints 
for the SED fitting, since the RW06 models may pre- 
dict variations of more than an order of magnitude in 
24 /jm flux density even when the IRAC flux densities are 
well-measured. Detections and upper limits were defined 
as follows: The formal uncertainty a on the aperture 
photometry is generally dominated by the uncertainty 
on the photometric background determination but also 
takes into account photon counting statistics (including 
the contribution of the zodiacal background light). Tak- 
ing F' to be the extracted flux density, sources for which 
F 24 < 5cr were considered undetected at 24 /zm, and up- 
per limits of F 24 < 5cr + F 24 for F 24 > and F 24 < 5er 
for ^4 ^ were assigned. The results of our aperture 
photometry are presented in Table 2, in which the 24 /zm 
quality flags (column 20) are set to 1 for detections and 
3 for upper limits. No photometric information can be 
obtained for sources falling within regions of the 24 /im 
mosaic where the nebular emission saturates (for exam- 
ple, in the vicinity of r\ Car); these are given quality 
flags of 0. For the 24 /zm detections we set the mini- 
mum uncertainty to 0.15F 2 4 to account for systematic 
uncertainties. One potential source of systematics is the 
mid-IR extinction law for A > 10 /im, which appears to 
vary with environment. Recent observational work with 
Spitzer has found that, in molecular clouds, the absorp- 
tion at 24 /im is similar to, not less than, the absorption 
at 8 /im (e.g. Flaherty et al. 2007). We therefore used 
the high-column-density (Ak > 1) mid-IR extinction law 
of McClure (2009) when fitting the RW06 models to the 
SEDs of the candidate YSOs. 

2.4. MSX Detections of Luminous Candidate YSOs 

We incorporate 8-21 /im photometry from the MSX 
Galactic Plane Survey Point Source Catalog (Price et al. 
2001) into our YSO sample in 2 ways: (1) Spatial cor- 
relation of MSX sources with the candidate YSOs and 



(2) visual identification of bright mid-IR sources within 
the CCCP field that were excluded from the Vela-Carina 
Catalog because they were either marginally resolved or 
saturated. Moderately saturated sources were included 
in the more complete Vela-Carina Point Source Archive, 
and we added 3 MSX sources with Archive matches to 
the sample. Only 11 PCYC sources have MSX counter- 
parts (indicated in Table 2), but these include the most 
luminous YSOs in the cloud. Most of these objects were 
previously cataloged by Rathborne et al. (2004) or Smith 
& Brooks (2007). The addition of up to 4 MSX data- 
points provides strong constraints on the mid-IR SEDs 
along with replacements for Spitzer photometry suffering 
from saturation at 8.0 or 24 /im. 

2.5. Removal of Contaminating IR Excess Sources 

While it is dominated by YSOs, the sample of reli- 
able IR excess sources (Table 1) may also contain the 
following types of contaminants: variable stars, dusty 
asymptotic giant branch (AGB) stars, unresolved plane- 
tary nebulae, and background galaxies (Robitaille et al. 
2008; Povich et al. 2009). Although few in number, lumi- 
nous AGB stars are the most important contaminants, 
as they can masquerade as massive YSOs. The majority 
of AGB stars have [8.0] - [24] < 2.2 mag, and the few 
extreme AGB stars that are redder are highly luminous 
and would be conspicuously bright (Whitney et al. 2008; 
Povich ct al. 2009). Applying this color cut removed 1% 
of sources. We correct for the remaining contamination 
by analyzing the spatial distribution of IR excess sources. 

While YSOs associated with the Carina Nebula ex- 
hibit a high degree of clustering (SP10), contaminat- 
ing sources, including foreground YSOs, are distributed 
uniformly. We used the nearest-neighbor group-finding 
algorithm described by Povich et al. (2009) to identify 
sources exhibiting a significant degree of clustering with 
respect to 3 "control" fields located outside of the CCCP 
field. The algorithm was tuned to be sensitive to group- 
ings of >10 sources. The control fields were 0?5 x 0?5 
boxes centered at (I, b) = (286?4, -0?9), (286?6, 0?2), and 
(288?3,0?1), containing a mean surface density of reli- 
able IR excess sources of S C on = 200 deg~ 2 . The surface 
density of reliable IR excess sources inside the 1.42 deg 2 
CCCP field is S CCCP = 1150 deg" 2 (Table 1), hence 
the expected contamination fraction from sources unas- 
sociated with the Carina Nebula is S con /Scccp = 17%. 
After removing sources from the sample that did not sat- 
isfy the clustering criterion, the source density in the final 
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PCYC is X£?£,q P — 1000 deg 2 , indicating a residual con- 
tamination level of Sq " cp — (E CCCP — S con ) = 50 deg -2 , 
or 5%. 

3. MODEL-BASED YSO CHARACTERIZATION AND 
CLASSIFICATION 

Basic IR photometric data for the 1439 YSOs in the 
PCYC are presented in Table 2. Using the Robitaillc 
et al. (2007) SED fitting tool, we fit YSO models from 
RW06 to the available photometry of each source. We 
allowed the fitting tool to accept a range of distances of 
2.3±0.05 kpc, guided by the well-established distance to 
77 Car from expansion measurements of the Homunculus 
nebula (Smith 2006b) and by the assumption that the 
depth of the Carina Nebula along the line of sight is not 
greater than its ~50 pc diameter observed on the sky. We 
define the set i of well-fit models for each YSO according 
t° Xi — Xo < 2iVdata, and assign a x 2 -weighted, normal- 
ized probability Pj to each model (Povich et al. 2009). 
This enables us to construct probability distributions 
for the model parameters and thereby constrain phys- 
ical properties such as luminosity, mass, and accretion 
rate. In general, the parameters are not normally dis- 
tributed, so the physical properties for individual sources 
must be interpreted with care. We present results from 
the model fitting to the PCYC sources in columns (2)-(8) 
of Table 3. For each YSO, we list characteristic values 
of 3 parameters along with corresponding uncertainties. 
The characteristic values for the stellar parameters mass 
(X = M+) and bolometric luminosity (X — Z/bol) ar e 
calculated as the probability-weighted means of the pa- 
rameter distributions, 

Circumstellar parameters in the RW06 models, such as 
envelope accretion rate (Y = M cnv ), span many orders 
of magnitude. We define the characteristic circumstellar 
parameters as the medians of the parameter distribu- 
tions, 

oo = Ml/2 (n (2) 

because this approach is less affected by any extreme 
outlying values returned by the fitting procedure. For 
both stellar and circumstellar parameters, we compute 
the la uncertainty <j(X, Y) as the standard deviation on 
the characteristic parameter value (X, Y) . 

The YSO models are divided into evolutionary stages 
that parallel the well-known empirical T Tauri classifica- 
tion system based on IR colors or spectral indices: Stage 
0/1 YSOs are still embedded in their infalling, natal en- 
velopes; Stage II YSOs have optically thick circumstellar 
disks; Stage III YSOs have optically thin disks (RW06). 
We construct probability distributions of evolutionary 
stage for each YSO and classify a YSO as Stage 0/1, 
II, or III if ^Pi(Stage) > 0.67; we refer to the classifi- 
cation of sources not meeting this criterion as "Ambigu- 
ous" (Povich et al. 2009). We prefer the RW06 stage 
classification system to the traditional T Tauri class sys- 
tem because (1) it makes the model-dependent nature of 
the taxonomy explicit rather than implicit, (2) the incor- 
poration of Ambiguous classifications acknowledges the 
limitations of imposing distinct categories on a continu- 
ous evolutionary sequence, and (3) a model-based system 



provides physical insight into intermediate- and high- 
mass YSOs for which a classification based on the empir- 
ical colors of low-mass T Tauri stars may not be mean- 
ingful. The distribution of stage classifications across the 
PCYC is presented in Table 1, and the classification for 
each YSO is given in Table 3 (column 8). 

The 410 YSOs detected in X-rays are identified in Ta- 
ble 3 by CCCP catalog name (column 9). Feigelson et 
al. (2011, hereafter Fll) divided the CCCP field into 
4 broad regions (A, B, C, and D) and also identified 
20 X-ray clusters and 31 smaller groups. We assign 
every PCYC source to one of these regions and possi- 
bly also to one of the clusters/groups according to the 
rcgion/cluster/group membership of the nearest CCCP 
catalog source (columns 10 and 11 of Table 3). 

The PCYC (Tables 2 and 3) is a compilation of in- 
dividual candidate YSOs in Carina that we hope will 
provide a rich dataset for future follow-up studies. For 
the present work, we are concerned primarily with the 
ensemble properties of the population. The large sam- 
ple size provides robust statistics to compensate for po- 
tentially large uncertainties on individual YSOs whose 
physical properties, such as stellar mass and bolometric 
luminosity, may not be well-constrained by the RW06 
models. 

4. GLOBAL PROPERTIES OF THE CARINA YSO 
POPULATION 

4.1. Spatial Distribution 

The positions of sources in the YSO catalog are plotted 
on an IRAC 8.0 fim mosaic image from the Vela-Carina 
survey in Figure 1. While YSOs are distributed widely 
throughout the CCCP field, the population generally 
traces the broad North-South stellar "backbone" join- 
ing the major ionizing clusters Trumpler (Tr) 14 and 16 
with the Collinder (Cr) 228/South Pillars region (Towns- 
ley et al. 2011, Fll). Contours of the major molecu- 
lar clouds traced in CO emission and the locations of 
8 dense C 18 molecular cloud cores (Yonekura et al. 
2005) are also plotted in Figure 1. While the spatial 
distribution of YSOs throughout the Carina Nebula ex- 
hibits highly complex structure, with clustering on multi- 
ple spatial scales, a general pattern is readily discernible. 
The majority of YSOs are located inside H II region cav- 
ities near, but less frequently within, the boundaries of 
dense molecular clouds and the ends of pillars. 

Another view of the PCYC is presented in Figure 2, in 
which a Vela-Carina 3.6 /im mosaic is presented with 
YSOs color-coded by mass. YSOs in 3 mass ranges 
roughly corresponding to the mass ranges of Herbig Be 
stars, Herbig Ae stars, and T Tauri stars are colored 
red, yellow, and green, respectively. Density contours of 
CCCP sources from Fll depict all of the major clus- 
ters as well as numerous smaller stellar groupings. The 
largest concentration of intermediate-mass YSOs is found 
not in the South Pillars, but in Tr 14 itself. Numerous 
sub-clusters of Tr 16 (Fll) also have associated YSOs. 

Several caveats regarding the Spitzer point-source de- 
tections must be kept in mind when analyzing the PCYC. 
The strong and variable nebular emission across the Ca- 
rina Nebula, especially at 8.0 and 24 ^m, creates com- 
plicated, position-dependent variations in the sensitiv- 
ity limit for mid-IR point-source detections. In addition, 
the IRAC images are confusion-limited in the dense clus- 
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Figure 1. Spitzer/IRAC 8.0 /im mosaic image with positions of the 1439 YSOs in the Pan-Carina YSO Catalog (PCYC) overplottcd. 
Colors identify YSO evolutionary stages: red=Stage 0/1, yellow=Stage II, green=Stage III, cyan=ambiguous. YSOs with X-ray counterparts 
are marked by crosses. Contours are CO integrated intensity over the velocity range Vlsr = — 30 km s — 1 to —10 km s — 1 , and the large 
magenta circles mark C ls O molecular cores (Yonekura et al. 2005). The solid green outline is the boundary of the CCCP survey area, 
while the dashed green lines enclose the area where the South Pillars observation of SP10 overlaps the CCCP. The bright, saturated star 
near the image center is r\ Car. Five well-known massive star clusters discussed in the text are labeled. 



ter centers. The PCYC is very incomplete near 77 Car, 
a highly saturated Spitzer source visible at the center 
of Figures 1 and 2, and in the Treasure Chest, a very 
young, embedded cluster (Smith et al. 2005), where only 
3 candidate YSOs are detected due to extreme crowd- 
ing and nebulosity. Other potential areas of sensitivity 
problems include the crowded center of Tr 14 and the 
ionization front between Tr 14 and the Car I molecular 
cloud core to the west (Yonekura et al. 2005; Ascenso 
et al. 2007) where the diffuse mid-IR nebular emission is 
bright. Nevertheless, the dearth of active star formation 
within the Car I cloud core appears to be real, not simply 
a sensitivity effect. 

All of the Spitzer YSO sub-clusters cataloged by SP10 
are apparent in Figures 1 and 2. In addition, we de- 
tect 3 new, tight clusters or groups of ~10 YSOs in the 



Northwest corner of the CCCP area, 2 of which are as- 
sociated with significant overdensities of X-ray sources 
(Fll). We discuss these new YSO clusters in more de- 
tail in §6. YSOs are also detected embedded in the ends 
of most dust pillars, both in the South Pillars and else- 
where in the Carina Nebula. Some of these are the driv- 
ing sources of the Herbig-Haro jets found in Hubble Space 
Telescope images (Smith et al. 2010a). YSOs are also as- 
sociated with several of the candidate proplyds identified 
by Smith et al. (2003), specifically proplyds 104405.4- 
592940 (PCYC 429), 104619.7-595044 (1126), 104632.9- 
600353 (1173), and 104519.3-594423 (842 and 841). All 
of the YSOs associated with proplyds have either Stage 
0/1 or Ambiguous classifications, and none are detected 
in X-rays. 
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Figure 2. Spitzer-IRAC 3.6 (im mosaic image with positions of the 1439 YSOs in the PCYC overplotted. Colors identify YSOs by mass 
range: red, (M*) > 3.1 Mq; yellow, 2 M©< (M*) < 3.1 Mq; green, (M*) < 2 Mq. Contours are from Fll and show the surface density 
of X-ray sources. The green outlines and cluster labels are the same as in Figure 1. 



4.2. YSO Mass Function and Present-Day Star 
Formation Rate 

Summing the individual probability distributions of 
M* defined by the set of models fit to each YSO, we 
construct the YSO mass function (YMF; Shepherd et al. 
2007; Povich et al. 2009; Kang et al. 2009, PW10) for all 
1439 YSOs in the PCYC (Figure 3). The YMF is de- 
fined similarly to the stellar initial mass function (IMF) , 
$(logm) = dN/d\ogm, where dN is the number of stars 
in the (logarithmic) mass interval (log m, log m+d log m) . 
The Carina YMF exhibits a power-law form ^(m) oc 
m -r Y MF f or ]\j c < m < 10 Mm, departing from the 



power-law due to incompleteness for m < M p « 3.1 M c 
(Figure 3a). The power-law slope Tymf = 3.2 ± 0.3 is 
significantly steeper than the standard Salpeter-Kroupa 
IMF slope Timf = 1-3 (Kroupa 2001). PW10 found 
a very similar YMF shape, with Tymf = 3.5 ± 0.6, 
for a sample of 488 YSOs in M17 SWex, a CMC com- 



plex extending ^50 pc outward from the bright Galactic 
H II region M17. 

In Figure 3b, the YMF has been broken down by most 
probable evolutionary stage. This shows that the overall 
shape of the YMF is dominated by a narrow component 
from Stage II objects, but also includes a broader com- 
ponent from Stage 0/1 objects. The YMF shape of the 
Ambiguous sources is a blend of these two components. 
(Stage III YSOs constitute only 6% of the PCYC and are 
not plotted.) Ambiguous YSOs dominate the (incom- 
plete) low-mass range, which is expected because fainter 
sources are more likely to lack 24 /jm detections and other 
photometric measurements of sufficient quality to tightly 
constrain their properties. Stage 0/1 YSOs are preferen- 
tially detected over Stage II YSOs at lower masses be- 
cause the reddest objects are generally brighter in the 
mid-IR for a given mass. The steep intermediate-mass 
slope of the YMF traces both Stage II YSOs and Am- 
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Figure 3. (a) YMF plot for all 1439 YSOs in the Catalog (his- 
togram points with error bars). Two IMFs are overplotted: power- 
law fit to the YMF (red) over the range m > M p and Kroupa 

(2001) IMF scaled to match the YMF at M p = 3.1 M Q (dashed 
red line). The YMF from PW10 for 488 YSOs in the M17 SWex 
cloud is overplotted for comparison (gray histogram), and its crit- 
ical mass M c = 3.9 Mq (see text) is indicated (dashed gray line). 
(6) Carina YMF binned up by a factor of 2 and subdivided by 
evolutionary stage. 

biguous YSOs that are most likely dominated by Stage II 
objects in this mass range, given the overall distributions 
of Stage 0/1 versus Stage II YSOs. This supports the idea 
that disk evolution is an important selection effect that 
steepens the observed YMF slope (PW10). 

SP10 considered the idea that rapid disk destruction in 
massive stars could preferentially remove them from the 
YSO sample and affect the observed mass function but 
concluded that this effect does not significantly impact 
an IMF derived from the YSOs in Carina. This con- 
clusion was based upon a perceived similarity between 
the empirical K$ and [3.6] luminosity functions of the 
Carina YSOs and the scaled K- and L-band luminosity 
functions (KLFs and LLFs) for the Orion Nebula Clus- 
ter (ONC; Muench et al. 2002). However, Muench et al. 
(2000) caution that LFs are not reliable tracers of IMFs 
when the photometry is dominated by IR excess emis- 
sion. To illustrate the problem, consider the bright end 
of the KLF. The few most luminous YSOs in the SP10 
sample, which are also included in the PCYC, have simi- 
lar Mk to 8 1 Ori C, the most luminous star in the ONC. 
But while the most massive YSOs in Carina are ~10 M Q 
objects, 9 1 Ori C is an O star of -40 M ! While the 
ONC is a young cluster, the majority of its stars have 
lost their inner disks, and the KLF used by Muench et al. 

(2002) to derive the cluster IMF was dominated by emis- 
sion from stellar photospheres. To better understand the 



IMF traced by the YSOs in Carina, we rely instead on 
the RW06 models to relate statistically the observed IR 
emission to stellar mass. Our approach is subject to its 
own set of biases and systematics, but it has been de- 
signed for the problem at hand and has been tested in 
numerous other studies (e.g. Whitney et al. 2008; Povich 
et al. 2009; Kang et al. 2009). We are also aided by the 
ability to compare Carina to M17 SWex, where a similar 
sample of YSOs has been analyzed using the same tools 
(PW10). 

The YMF of 488 YSOs in M17 SWex from PW10 is 
overplotted in Figure 3a for comparison with the PCYC. 
The similarity between these YMFs is striking; the only 
obvious difference appears to be a lower mass complete- 
ness limit for the PCYC. Both Carina and MI7 are 
100-pc scale CMC complexes at similar heliocentric dis- 
tances, each containing several 10 5 M Q in molecular gas 
mass, and each in the process of forming a large OB asso- 
ciation. The sensitivity of the GLIMPSE observations to 
YSOs in M17 SWex was limited primarily by interstellar 
extinction through the IR dark cloud. In contrast, the 
interstellar extinction to the Carina YSOs is much lower, 
and the sensitivity is limited instead by bright nebular 
emission produced by the H II regions. It is not obvious, 
however, why high extinction should be a much greater 
impediment to YSO detections than bright nebulosity. 
The implied difference in sensitivity is dramatic, mean- 
ing that if the YMF of PW10 is intrinsically similar to the 
PCYC, its completeness at m — 3.1 M (corresponding 
to main-sequence B9 stars) is as low as ~40%. 

The M17 SWex CMC complex appears to be form- 
ing its first generation of stars and lacks very massive, 
early O stars (PW10). The Carina molecular cloud com- 
plex has been thoroughly disrupted by the process of 
forming tens of thousands of stars, including some of 
the most massive stars known in the Galaxy (Smith & 
Brooks 2007). The environmental differences between 
these two massive star-formation regions lead us to ex- 
pect that there are intrinsic differences between the YSO 
populations, and these differences might be manifest in 
the YMFs. In particular, we expect that the PCYC in- 
cludes a larger proportion of more evolved YSOs than 
the MI7 SWex population. 

A useful empirical comparison of the PCYC with the 
M17 SWex YSOs is provided by the plot of the IRAC 
[3.6] — [4.5] versus [4.5] — [5.8] color space presented in 
Figure 4. PW10 identified a class of predominantly Stage 
0/1 YSOs that exhibited excess 4.5 /im emission (denoted 
[4.5]E) with respect to the RW06 model SEDs and hy- 
pothesized that such sources are unresolved analogs of 
the extended green objects (EGOs), candidate molec- 
ular outflows from massive YSOs cataloged from the 
GLIMPSE images by Cyganowski et al. (2008). In this 
interpretation, such sources are expected to be very 
young, with ages <10 5 yr, since outflows are powered by 
rapid accretion (Richer et al. 2000; Shepherd et al. 2007). 
We find that -50% of [4.5]E candidates from PW10 arc 
found above the line 

[3.6] - [4.5] > l.I9([4.5] - [5.8]) + 0.5 (3) 

in Figure 4; these "strong" [4.5]E sources represent 7% 
of the YSOs in M17 SWex, but <1% of the PCYC. SP10 
searched the Spitzer images of Carina for EGOs and 
found only four candidates. The corresponding dearth of 
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Figure 4. IRAC color space enabling the separation of YSOs with 
excess 4.5 fim ([4.5]E) emission. PCYC sources detected in these 
3 bands are plotted along with sources from the M17 SWex YSO 
sample (PW10), with [4.5]E YSOs plotted in green. Candidate 
[4.5]E objects are located near or above the dashed line. Reddening 
vectors for Ay = 30 mag using the extinction laws of Weingartner 
& Draine (2001) and Indebetouw et al. (2005) are shown as open 
and filled arrows, respectively. The cross represents the typical 
photometric uncertainty. 

[4.5]E sources supports the conclusion that these objects 
trace the same phenomenon as the EGOs. The differ- 
ence between the two regions could be partly environ- 
mental; SP10 suggest that the harsh radiation fields per- 
meating the Carina Nebula may dissociate molecules in 
the outflows, removing the source of the EGO emission. 
But this comparison also indicates that very young YSOs 
with high accretion rates are less frequent in Carina than 
in M17 SWex. Indeed, only 17% of the PCYC are un- 
ambiguous Stage 0/1 YSOs, compared to 27% Stage 0/1 
YSOs in the M17 SWex sample (PW10). 

The presence of a dusty, circumstellar disk (or an enve- 
lope that will presumably evolve into a disk) is a prereq- 
uisite for inclusion in the PCYC, hence the upper limit 
on the age of the YSOs is approximately the disk life- 
time. The "canonical" disk lifetime for a solar-mass, 
classical T Tauri star is ~2 Myr (Haisch et al. 2001). If 
the primary mechanism for destroying dust disks around 
intermediate-mass YSOs is photo-evaporative dissipation 
driven by radiation from the central star (Hollenbach et 
al. 1994; Monnier & Millan-Gabet 2002), then the disk 
lifetime should decrease as a strong function of stellar 
mass. Disk lifetimes for intermediate-mass stars, includ- 
ing Herbig Ae/Be stars, appear to be significantly shorter 
than for T Tauri stars, typically <1 Myr (Hernandez 
et al. 2007; Povich et al. 2009, PW10). The disk frac- 
tion among intermediate-mass stars may therefore be one 
of the most sensitive chronometers available for massive 
star-formation regions like Carina. 

PW10 showed that the shape of the YMF depends crit- 
ically on the mass-dependent disk fraction, and hence the 
recent star formation history. Assuming a constant SFR, 
the observed YMF slope is steepened by disk evolution 
above a critical mass m > M c for which the disk lifetime 
Td{m) is less than the age spread of the YSO population 
To. For m < M c , Td{m) > tq, the disk fraction becomes 
unity, and the YMF traces the underlying IMF, presum- 
ably a Salpeter-Kroupa slope (Fymf — > Timf = 1.3). 
As a result, M c decreases with time. Unfortunately, for 



the PCYC M c < M p = 3.1 M Q , meaning that the crit- 
ical mass falls in the range where our observations are 
incomplete. This limit means that the critical mass of 
the PCYC is lower than M c = 3.9 M Q found by PW10 
for the M17 SWex sample (Figure 3a), and this could 
be explained by disk evolution proceeding over a longer 
timescale in the Carina complex. The PCYC does con- 
tain a significant population of Stage II/III YSOs with 
relatively evolved disks. This population is located near 
the line [3.6] — [4.5] = in Figure 4, hence reddening by 
interstellar dust is negligible, and the lack of IR excess 
emission at wavelengths <4.5 fim suggests that the inner 
disks have been cleared of dust. 

Assuming that the underlying IMF of the Carina young 
stellar population has a normal Salpeter-Kroupa form, 
we can estimate the present-day SFR in Carina from the 
PCYC. Scaling the IMF of the ONC (Muench et al. 2002) 
to match the PCYC YMF where it departs from a power- 
law at M p — 3.1 Mq gives a prediction of >2 x 10 4 
YSOs to a limiting mass of 0.1 M Q , corresponding to a 
total mass in stars of >1.6 x 10 4 Mq. These numbers 
are lower limits due to incompleteness. Adopting a con- 
servative upper limit of r < 2 Myr as the age spread 
sampled by the PCYC, we arrive at a lower limit on the 
recent SFR in the Carina Nebula of >0.008 M Q yr" 1 . If 
the Galactic SFR is a few M Q yr _1 (Murray & Rahman 
2010; Robitaille & Whitney 2010), then the star forma- 
tion activity of the entire Milky Way is equivalent to a 
few hundred Carina complexes. 

5. PROPERTIES OF THE X-RAY-EMITTING YSO 
POPULATION 

The PCYC contains 410 X-ray-detected, predomi- 
nantly intermediate-mass YSOs, including 62 Stage 0/1 
YSOs that are candidate X-ray-emitting protostars, the 
largest such sample compiled to date (Table 1). Because 
Herbig Ae/Be stars fall into an X-ray "desert" between 
low-mass T Tauri stars with convection-driven emission 
and OB stars with wind-driven emission, we expect that 
many of the PCYC sources are intrinsically X-ray quiet. 
This motivates us to assume the null hypothesis that 
X-ray emission apparently associated with intermediate- 
mass YSOs is generally produced by an unresolved low- 
mass, T Tauri binary companion (Stelzer et al. 2006; 
Evans et al. 2011). In such a scenario, the intermediate- 
mass primary is X-ray quiet but dominates the IR emis- 
sion, and the physical properties derived from the IR 
SED modeling are expected to be uncorrelated with the 
presence of an X-ray-detected secondary. This null hy- 
pothesis can therefore be disproved if we observe statisti- 
cally significant differences between the modeled proper- 
ties of YSOs with X-ray counterparts and those without. 

The detection of a faint X-ray source is highly depen- 
dent on location, because the sensitivity of the ACIS 
observations varies significantly across the CCCP field 
(Broos et al. 2011a). The sensitivity limit of the PCYC 
has a complicated positional dependence due to large 
spatial variations in extinction and nebulosity. Physi- 
cal properties inferred from the IR modeling can range 
over orders of magnitude, and observed X-ray luminos- 
ity Lx versus stellar mass (or Lx versus I/boi) relations 
have over an order of magnitude scatter (e.g. Preibisch 
et al. 2005). Potential correlations between properties 
derived independently from X-ray and IR observations 
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Figure 5. Same as Figure 3, except only the 410 YSOs with X- 
ray counterparts are included in the X-YMF. The gray histogram 
shows the full PCYC YMF from Figure 3. 

are likely to be overwhelmed by severe noise. To avoid 
this problem, we divide the PCYC into CCCP-matched 
(PCYC-X) and unmatched (PCYC-U) subsamples and 
compare physical properties derived from SED fitting 
between these subsamples. The PCYC-U includes both 
YSOs that are intrinsically X-ray quiet and YSOs that 
happened to fall below the CCCP detection limit. There 
is likely to be substantial overlap in the physical pa- 
rameters between these subsamples. Rather than at- 
tempt the potentially intractable task of correcting for 
the incredibly complicated, competing sensitivity varia- 
tions between the X-ray and IR sample selections, we 
regard sample overlap as an inevitable, additional source 
of noise in our comparisons. At best, we might hope to 
discern some general, qualitative trends from the ensem- 
ble population, given the large sample size. 

The YMF of the 410 YSOs in the PCYC-X (the X- 
YMF) is plotted in Figure 5. The intermediate-mass 
slope Tymf = 3.2 ± 0.6 and power-law departure mass 
M p = 3.1 M Q (Figure 5a) are consistent with the full 
PCYC YMF. This similarity is striking, given that we 
might have expected the X-YMF to skew toward lower- 
mass YSOs drawn from the T Tauri population. Instead, 
there appears to be a subtle trend in the opposite direc- 
tion, with the m < 2 Mq mass range suppressed in the 
X-YMF. This trend is also apparent when the X-YMF is 
plotted by evolutionary stage (Figure 5b). 

Differences between the PCYC-X and PCYC-U are 
easier to discern from the cumulative mass distributions 
(Figure 6), constructed from the characteristic (M*) for 
each YSO (column 4 of Table 3). Two curves are plotted 
in each panel, representing the PCYC-U and the PCYC- 
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Figure 6. Cumulative mass functions comparing X-ray detected 
(red) and undetected (black) YSOs: (a) All YSOs in the PCYC, 
(b) Stage II YSOs, and (c) Stage 0/1 YSOs. The number of YSOs 
in each curve is printed in each panel, and each panel header gives 
the probability from a KS test that the two samples are drawn 
from the same underlying distribution. The vertical dash-dotted 
lines indicate the YMF power-law departure mass, M p = 3.1 Mq 
(see Figure 3). 



X (red) subsamples. The panel headers give the prob- 
ability P(KS) from a Kolmogorov-Smirnov (KS) test 
that the two curves are drawn from the same underly- 
ing distribution. Figure 6 a shows that the PCYC-X as 
a whole includes a significantly lower fraction of sources 
with m < 2 Mq. Figure 6b reveals that this difference 
is produced by the Stage II sources. Evidently the Vela- 
Carina Spitzer observations are more sensitive than the 
CCCP observations to low-mass, disk-bearing stars. A 
KS test shows no significant difference between the cu- 
mulative mass distributions of the two samples of Stage 
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0/1 sources, although the two curves do appear to diverge 
for m > M p (Figure 6 c). 

5.1. Intermediate-Mass, Disk- Dominated YSOs with 
X-ray Counterparts 

In Figure 66 the mass distributions of both the PCYC- 
X and PCYC-U samples appear to be very similar for 
m > 2 M©. In Figure 7 we present a comparison of 6 
key YSO model parameters for the unambiguous Stage II 
YSOs in the intermediate-mass range 1.6 M Q < (M*) < 
4 M©, equivalent to A through late B stars on the main 
sequence. The parameters plotted are stellar mass (M*), 
bolometric luminosity (L bo {) (which includes luminosity 
derived from accretion), stellar temperature (7*), nor- 
malized envelope accretion rate (M env /M*), disk mass 
(Mdisk), and total extinction to the stellar surface Ayj 
(defined as the sum of the contributions to extinction 
produced by dust along the line of sight Ay and circum- 
stellar dust A v _ c ). The characteristic values (X) of the 
stellar parameters X = (M*, L bo \, T*) were calculated 
using Equation (1), and the characteristic values of the 
circumstellar/environmcntal parameters (Y) were calcu- 
lated using Equation (2). 

We further subdivide our intermediate-mass Stage II 
sample by stellar temperature, into sets with (T*) < 
5500 K (167 YSOs, 77 in the PCYC-X) and (T*) > 
5500 K (195 YSOs, 83 in the PCYC-X). We were mo- 
tivated to make this division by the hypothesis that 
intermediate-mass stars with cooler photospheres are 
more likely to have convective atmospheres and per- 
haps generate X-rays by a mechanism analogous to that 
of lower-mass T Tauri stars (e.g. Stassun ct al. 2004; 
Preibisch et al. 2005). 

We find that none of the mass, luminosity, or temper- 
ature distributions (Figures 7abc) shows any significant 
correlations with the presence of X-ray emission. The 
normalized envelope accretion rates are low by definition 
for Stage II YSOs, and the distributions shown in Fig- 
ure 7 confirm that these objects lack significant infalling 
envelopes. We do observe that the distribution functions 
for the low-T* division (dashed curves) have systemati- 
cally lower luminosity and mass compared to the high-T* 
division, which is expected from PMS stellar evolution- 
ary tracks. 

There are no statistically significant differences in the 
cumulative distribution functions of the model parame- 
ters among the high-T* division that correlate with X- 
ray emission. 17 This result is consistent with our null 
hypothesis; we cannot rule out the possibility that X- 
ray emission associated with Stage II YSOs with (T*) > 
5500 K actually originates in unresolved, lower-mass T 
Tauri companions. In contrast, the distribution func- 
tions of Mdisk and Ay,t (Figures 7ef) definitively fail the 
KS test for the low-T* division and demonstrate that X- 
ray emission is not a completely random phenomenon 
among the intermediate-mass Stage II YSOs. PCYC-X 
sources are significantly underrepresented among YSOs 
with the lowest disk masses, (Mdisk) < 10~ 3 M© (Fig- 

17 The one exception is normalized envelope accretion rate, 
P(KS) = 0.00 in Figure 7d\ however this is not a meaningful re- 
sult because the envelope accretion rates are near zero in Stage II 
YSOs, and this parameter is not expected to be well constrained 
by the SED fitting. 



ure 7e) . We have uncovered a population of intermediate- 
mass PMS stars with massive disks and intrinsic X-ray 
emission, but the X-ray emission apparently dies off as 
the disks evolve. This result seems diametrically opposed 
to the observational paradigm in which X-ray luminos- 
ity from low-mass PMS stars increases as these objects 
evolve from classical T Tauri stars to weak-lined T Tauri 
stars (Preibisch et al. 2005; Tellcschi et al. 2007), a tran- 
sition that reflects the cessation of disk-fueled accretion 
that may suppress X-ray emission (Stassun et al. 2004). 

The main result illustrated by Figure 7 is that disk 
mass correlates far more strongly with X-ray emission 
than any of the stellar parameters, demonstrating that 
rapid disk evolution provides a highly sensitive chronome- 
ter for intermediate-mass PMS stars. The apparent 
paradox with T Tauri X-ray emission is resolved when 
we recall that this trend is observed only among our 
carefully selected sub-division of intermediate-mass PMS 
stars with (M*) > 1.6 M© and (T*) < 5,500 K. These 
stars will presumably evolve into Hcrbig Ac/Be stars, at 
which point their X-ray brightness would fall below the 
CCCP detection limits (Stelzer et al. 2006; Broos et al. 
2011a). At present, they are very young objects with rel- 
atively massive disks, still contracting toward the main 
sequence, and they have spectral types of G2 or later. 
They could thus be regarded as more massive, more lu- 
minous analogs of classical T Tauri stars. 

A direct imprint of circumstellar disks on X-ray emis- 
sion can be inferred from the distribution of total ex- 
tinction (Ay t ; Figure If). The top 5% of the PCYC-U 
reach absorptions up to ^100 times higher than the top 
5% of the PCYC-X. These sources are probably high- 
inclination systems, where nearly edge-on disks com- 
pletely absorb any X-ray emission along the sightline to 
the central stars. 

5.2. Candidate X-ray- Emitting Protostars 

The 62 unambiguous Stage 0/1 YSOs in the PCYC-X 
subsample are presented in Table 4, which lists values for 
some of the YSO model parameters along with 3 observed 
X-ray quantities; total-band (0.5-8 keV) net counts, me- 
dian energy (£ mo dian), and energy flux (F t ). As Fig- 
ure 6c shows, the mass distributions of Stage 0/1 YSOs 
are very similar between the PCYC-X and PCYC-U sub- 
samples for (M*) < M p . To investigate which properties 
other than mass may correlate with X-ray emission from 
protostars, in Figure 8 we present distribution functions 
comparing the same 6 YSO model parameters used to 
analyze the intermediate-mass Stage II population (Fig- 
ure 7), this time plotting the unambiguous Stage 0/1 
YSOs with (M*) < 2.5 M©. 

KS tests reveal that the mass distributions for these 
low-mass Stage 0/1 YSOs are statistically the same be- 
tween the PCYC-X and PCYC-U (Figure 8a), yet the 
luminosity distributions arc significantly different, with 
the X-ray emitters systematically preferring lower values 
of Lboi (Figure 86). Hence we may deduce that X-ray 
detected Stage 0/1 YSOs are less luminous for a given 
stellar mass. Furthermore, both the normalized enve- 
lope accretion rates M env /M+ and total absorption Ay_ t 
(columns 6 and 7 of Table 4) show parameter distri- 
butions skewed to significantly lower values for the X- 
ray emitters, while is marginally higher for the X- 
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ray emitters. As protostars contract, they become hot- 
ter and less luminous, their accretion rates slow, and 
outflow-carved, bipolar cavities grow in their envelopes, 
reducing the obscuration to the central star. Taken to- 
gether, the distributions shown in Figure 8 demonstrate 
that brighter X-ray emission correlates with more evolved 
Stage 0/1 YSOs. 

Previous identifications of candidate X-ray emitting 
Class 0/1 YSOs in nearby star formation regions have 
observed characteristically hard X-ray spectra, with 
^median ^ 2.5 keV (Getman et al. 2007; Prisinzano et al. 
2008). X-ray spectra observed from embedded protostars 
are expected to be hard because soft X-rays are preferen- 
tially absorbed both by the dense circumstcllar envelopes 
and the surrounding dense molecular cloud cores (Get- 
man et al. 2007). About half (26/62) of our sample of 
candidate X-ray emitting Stage 0/1 YSOs have hard X- 
ray counterparts in the CCCP catalog (-E me dian > 2 keV); 
these objects are high-probability candidate X-ray emit- 
ting protostars, grouped together in Table 4. The re- 
maining Stage 0/1 YSOs in Table 4 arc matched to 
soft X-ray sources, more consistent with expectations 
for Class II/III objects (Preibisch et al. 2005; Prisin- 
zano et al. 2008). YSOs in this latter group may be 
in transition between the envelope-dominated and disk- 
dominated phases. Many of them would not have been 
selected as Class 0/1 objects by the conservative mid-IR 
empirical criteria of Prisinzano et al. (2008), for example. 
The model-based Stage 0/1 classification, however, indi- 
cates that the YSOs in Table 4 are still very young (they 
possess envelopes); whether or not they are identified 
as "candidate protostars" is primarily a semantic choice. 
As expected, the IR SED modeling indicates that the 
soft X-ray group is less embedded than the hard X-ray 
group, with median Ay,t = 10 mag compared to 30 mag 
(column 7 of Table 4). Recalling that most YSOs in 
the Carina Nebula are found outside of dense molecu- 
lar clouds (§4.1), the typical absorption may also be less 
than expected from observations of YSOs in other star- 
forming molecular clouds that have not been disrupted 
by feedback from extremely massive stars. 

6. BUTTERFLY COLLECTION: INTRIGUING YSO 
SUB-CLUSTERS 

In this section, we identify and briefly discuss YSOs 
found in candidate embedded sub-clusters in the CCCP 
field. These YSOs are highlighted, with their group 
membership identified, in Tables 2 and 3. There are a 
few caveats to remember. The properties of individual 
YSOs may not be well-constrained by the RW06 mod- 
els. Sources might be multiple or otherwise confused 
in our multiwavelength datasets, and apparent spatial 
matches between CCCP sources and Vela-Carina sources 
do not guarantee physical association. We believe that 
these clusters are sufficiently interesting individual "but- 
terflies" among the larger stellar population in the Carina 
complex to merit targeted follow-up multiwavelength ob- 
servations that would provide higher spatial resolution, 
wider wavelength coverage, and/or spectroscopy. 

6.1. A Perfectly Ordinary Stage 0/1 YSO Associated 
with an Extraordinary Obscured X-ray Cluster 

Townsley et al. (2011) report the discovery of a group 
of 4 X-ray sources in the South Pillars region that are 
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Figure 9. Top: SED of PCYC 699 (points) and well-fit YSO mod- 
els (solid curves). The heavy curve represents the best-fit model, 
and the dashed curve is the stellar photosphere of this model as it 
would appear in the absence of circumstellar extinction (removing 
A\/,c)- Bottom: Plot of parameters M* versus Av >c for the set 
of well- fit YSO models. The gray shaded area represents the full 
parameter space sampled by the RW06 YSO models. 

remarkably luminous, hard, and highly-obscured. Three 
of these sources were matched to a single Vela-Carina 
Catalog source, G287.9186-01.2913 (Broos et al. 2011a), 
or PCYC source 699 (Tables 2-4). The primary CCCP 
match to PCYC 699 is CXOGNC J104451.91-602511.9, 
the brightest (F t — 5 x 10 -13 erg cm~ 2 s _1 ) and second- 
hardest (-Emedian = 5.2 keV) X-ray source associated 
with any candidate protostar in the Carina Nebula (Ta- 
ble 4). This source is variable and spectacularly lu- 
minous in X-rays, with Lx = 0.5 L©, corrected for 
Ay,t = 250 mag total absorption (Townsley et al. 2011). 
In stark contrast, there is nothing remarkable about 
PCYC 699, save perhaps for the fact that its physical 
parameters are more poorly constrained than usual. The 
RW06 model fits to the SED of PCYC 699 arc plotted in 
Figure 9 along with a parameter plot showing versus 
Ay, c for each fit. A degeneracy in the long- wavelength 
SED beyond 24 fim, where we have no data, produces 2 
distinct families of model fits (Figure 9): intermediate- 
mass YSOs with relatively high circumstellar extinction 
and low-mass YSOs with lower circumstellar extinction. 
The model fits corresponding to the highest extinction 
are consistent with the absorption of the X-ray source, 
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but even taking into account the wide range of stellar 
mass and luminosity allowed by the model fits, there is no 
interpretation that would allow PCYC 699 to be consid- 
ered a massive YSO. It is a relatively faint 24 /xm source, 
and, unlike the most luminous YSOs in the PCYC, it is 
neither an MSX source nor an IRAS source. Thus we are 
left with a mystery surrounding the origin of the bright 
X-ray emission. The X-ray properties indicate a lumi- 
nous, obscured, compact massive star cluster, but there 
is no commensurate IR emission. 

6.2. Five Prominent YSO Sub-Clusters or Groups 

In Figure 10 we present Spitzer multicolor images 
zoomed in on 5 prominent YSO clusters, 3 of them new 
discoveries from the combined Spitzer /Vela-Carina and 
CCCP observations. Extending the nomenclature of the 
Spitzer clusters identified by SP10, we call the new clus- 
ters Spitzer P, Q, and R (Table 3). 

Spitzer P and Q: (a 2 ooo, ^2000) = (10' 1 43"\ -59°26'). 
These two clusters appear in projection to be close neigh- 
bors (Figure 10a), located ~8' northwest of Tr 14. Both 
contain a high fraction of unambiguous Stage 0/1 YSOs 
and show no sign of dynamical relaxation, suggesting 
extreme youth. The location of Spitzer P corresponds 
to a C 18 core (Figure 1), and its dense center con- 
tains a concentration of X-ray sources (Group 4; Fll). 
Spitzer Q was not designated an X-ray group by Fll, 
although it also coincides with an enhanced density of 
X-ray sources (Figure 2). Possibly the most massive X- 
ray-emitting Stage 0/1 YSO in the CCCP field is found in 
Spitzer Q. PCYC 179 is an extremely red Spitzer source, 
detected by MSX but not by 2MASS (Table 2), with 
(Lboi) ~ 4 x 10 3 L Q and a very high envelope accretion 
rate, (M onv ) - 2.5 x 10~ 3 M Q yr" 1 (Table 3). The as- 
sociated CCCP source is very hard, i? me dian = 4.0 keV 
(Table 4) . If the intrinsic X-ray emission is the soft spec- 
trum produced by a normal massive stellar wind (Lucy 
& White 1980), then it must be very luminous, with the 
soft counts almost completely absorbed. Alternatively, a 
more exotic emission mechanism may produce in intrin- 
sically hard X-ray spectrum, or the CCCP source is an 
unresolved, low-mass companion. 

Spitzer R: (a 20 oo, £2000) = (10 /l 44 m , -59°08')- In the 
extreme northern corner of the CCCP field we have dis- 
covered another very young, tight grouping of YSOs 
called Spitzer R (Figure 106). This cluster is roughly 
centered on a CO peak that was not detected as a C 18 
core by Yonekura et al. (2005, and see Figure 1), and it 
appears to be forming intermediate-mass stars but not 
massive stars. There is no associated X-ray cluster, but 
this may simply be a sensitivity effect given that Spitzer 
R lies close to the edge of the CCCP field and was ob- 
served far off-axis by ACIS, where the sensitivity is re- 
duced (Broos et al. 2011a). 

Cr 232: (a 2000 , <5 20 oo) = (10 ,l 44.5 m , -59°34'). Cr 232 
is a well-known young massive cluster to the west of Tr 
14, containing an 03.5 V((f)) star, HD 93250, and an 
05 V star, HD 303311 (Smith 2006a, and references 
therein). Cr 232 is a large X-ray cluster (Fll) appar- 
ently centered near the end of a broad dust pillar oriented 
toward r\ Car and Tr 16 (Figure 10 c). Cr 232 is un- 
doubtedly richer than the newly-discovered Spitzer YSO 
clusters. Its location suggests that its formation may 



have been triggered recently by feedback from the mas- 
sive stars in Tr 16. Apparently in the middle of the dust 
pillar associated with Cr 232 we detect a candidate X- 
ray-emitting, highly-embedded, massive Stage 0/1 YSO, 
PCYC 556 (Table 3), similar to PCYC 179 in Spitzer 
Q, but the associated X-ray source is fainter by 0.4 dex 
and has -Emedian = 1-9 keV, placing it among the soft 
X-ray group in Table 4. PCYC 556 is among the red- 
dest IRAC sources in the PCYC ([3.6] - [4.5] = 1.4 mag, 
[5.8] - [8.0] = 1.0 mag; Table 2). This relatively soft X- 
ray emission apparently originating from an embedded 
massive star suggests a highly luminous, intrinsically soft 
X-ray spectrum, perhaps powered by shocks in a strong 
stellar wind or accretion flow. Alternatively, the CCCP 
source may actually be a less obscured star that happens 
to be confused with the bright IRAC source. 

Spitzer O: (a 2000 , c5 2000 ) = (10 /l 45 m , -60°26.5')- Iden- 
tified recently by SP10, Spitzer O, located in the extreme 
southeast corner of the South Pillars, is a far outpost 
of massive star formation in the Carina complex (Fig- 
ure lOd). The cluster is associated with the southern 
of two C 18 cores in the molecular cloud known as the 
Giant Pillar (Figure 1). Although it is near the edge 
of the CCCP field, where the ACIS sensitivity is lower, 
Spitzer O is nevertheless associated with a surface den- 
sity enhancement in the X-ray sources (Figure 2), which 
was not listed by Fll on account of its remote location. 
Spitzer O includes a small IR dark cloud that hosts one 
of the few examples of an EGO found in the Carina Neb- 
ula (SP10). This EGO appears to be associated with 4 
Stage 0/1 YSOs (Table 3) in a tight, linear grouping. Ex- 
tending this line westward along the IR dark cloud (to 
the right in Figure lOd), we find a Stage II YSO and a 
candidate low-luminosity ultracompact H II region, sug- 
gesting a time-sequence of intermediate-mass star for- 
mation. Further to the west, Spitzer O also contains 
at least one massive star inside a compact H II region, 
identifiable as a small bubble outlined in 8 /im emission 
containing diffuse 24 /im emission from dust heated by 
the central ionizing source, most likely a late O-type star 
(Pll). Two similar but slightly less bright H II regions, 
complete with their own ionizing stars, are visible to the 
north and south of cluster O. All 3 candidate OB stars 
in Figure lOd are detected in X-rays, and none appears 
to have been previously cataloged (Pll). 

7. DISCUSSION 

7.1. A Recent Star Formation History for the Carina 

Nebula 

In §4.2, we derived a lower limit on the recent, global 
SFR in the Carina Nebula of >0.008 M yr _1 , averaged 
over the past To < 2 Myr, traced by the IR excess sources 
included in the PCYC. Feedback from massive stellar 
winds and radiation has been invoked by many authors 
as the driver of ongoing star formation (e.g., Smith et 
al. 2000; Smith & Brooks 2007, SP10). If triggering by 
massive stars is indeed the dominant mode of continuing 
star formation in GMCs that host energetic H II regions, 
we can assess its relative importance by comparing the 
present-day SFR to the historical SFR that produced the 
massive star clusters that ionize the Carina Nebula and 
provide the source of the feedback. 
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Figure 10. Spitzer color images (red = 24 /im, green = 8.0 fim, blue = 4.5 /im) highlighting the 5 tight groupings or clusters of YSOs 
discussed in the text. In all panels, YSOs are marked by circles color-coded by most probable evolutionary stage as in Fig. 1, and crosses 
mark the positions of CCCP X-ray catalog sources. Where an X-ray source has been matched to a YSO, the crosses have been colored to 
match the YSO classification. Dashed ovals have been drawn subjectively to enclose the YSOs that appear associated with each grouping. 
Coordinate axes are equatorial (J2000). In panel (c), two cataloged, X-ray-emitting early O stars in Cr 232 are marked by orange diamonds 
and crosses; these are HD 93250 (left) and HD 303311 (right). Regions at the bottom edge of this image appear turquoise due to saturation 
of the 24 fira nebular emission. In panel (d), three candidate X-ray detected OB stars from Pll are marked by blue diamonds and crosses, 
and a candidate ultracompact H II region is highlighted with a cyan diamond. 



Smith & Brooks (2007) estimated a total stellar pop- 
ulation of 5-8 x 10 4 stars down to the hydrogen-burning 
limit by extrapolating from the known O-type stars in 
Carina (Smith 2006a). The relative rarity of O stars 
requires an enormous extrapolation over the IMF to 
reach the low-mass stars that contain the bulk of the 
mass (Kroupa 2001). Measurements of the IMF among 
low-mass stellar populations (from the Galactic field 
and open clusters) and high-mass stellar populations 
(from very massive Galactic and extragalactic clusters) 
have generally found slopes consistent with the Salpeter- 
Kroupa value, but region-to-region variations do exist 
(Bastian et al. 2010). Because very massive clusters with 
a well-populated upper IMF are located at such large dis- 
tances that individual low-mass stars generally cannot be 
resolved, there has been no simultaneous, self-consistent 
measurement of the IMF from solar-mass stars to the 
most massive stars. Compounding the problem, theo- 
retical calibrations of the masses of O stars encounter a 
factor of ^2 discrepancy between "spectroscopic" masses 



derived from surface gravity and "evolutionary" masses 
expected from stellar structure (Martins et al. 2005; Wei- 
dner & Vink 2010). Finally, in the specific case of the 
Carina Nebula the completeness of the sample of known 
O stars (Smith 2006a) may be as low as 50% (Pll). 

The CCCP X-ray catalog itself offers an independent 
measure of the total stellar population. The X-ray lumi- 
nosity function (XLF) has proven to be a reliable tracer 
of the IMF (Getman et al. 2006) that is most sensitive 
to low- and intermediate-mass PMS stars. Fll scale the 
XLF from the Chandra Orion Ultradeep Project (COUP; 
Feigelson et al. 2005), for which the underlying IMF is 
well-characterized (Muench et al. 2002), to the XLFs 
from 4 distinct broad spatial regions within the CCCP 
and predict a total population of ~10 5 stars above the 
brown dwarf mass limit. Because X-ray detection effi- 
ciencies are lower for YSOs versus diskless PMS stars 
(Preibisch et al. 2005) and for heavily obscured versus 
lightly obscured stars (e.g. Figure 3 in Feigelson et al. 
2005), the simple XLF scaling analysis implicitly assumes 
that any differences in the global disk fractions or the 
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Figure 11. Comparison of global total-band (0.5-8 keV), 
absorption-corrected X-ray luminosity functions (XLFs) for disk- 
less candidate PMS stars in the CCCP to the XLF (gray histogram) 
of 839 lightly obscured cool stars from the COUP sample (Fcigclson 
et al. 2005). XLFs arc analyzed separately for the lightly-obscured 
(^median $ 2 keV; solid black histogram) and heavily-obscured 
(^median > 2 keV; dashed histogram) CCCP samples, and the 
scale factors used to match the COUP XLF to each CCCP XLF 
are labeled on the plot. Best-fit (negative) power-law slopes to 
each XLF are ovcrplotted as lines; these are Txlf = 0-9, 1.3, and 
1.0 for the COUP, lightly-obscured CCCP, and heavily-obscured 
CCCP XLFs, respectively. 



distributions of X-ray absorption between COUP and 
CCCP are negligible. 

In Figure 11 we present a more detailed analysis of 
the CCCP XLF. Two XLFs, separated according to me- 
dian energy into lightly-obscured (-E mo dian < 2 keV) 
and heavily-obscured (-^median > 2 keV) samples, were 
constructed using all CCCP sources identified as clus- 
ter members that were sufficiently bright to enable the 
derivation of absorption-corrected luminosities Lx (Get- 
man et al. 2010; Broos et al. 2011b). The 410 CCCP 
sources matched to YSOs in the PCYC were excluded, 
as were ~200 X-ray detected known or candidate OB 
stars (Pll). The resulting XLFs (Figure 11) are clean 
samples of diskless PMS stars that can be more directly 
compared to the 839 lightly obscured stars in the COUP 
XLF, the vast majority of which lack circumstellar disks 
(Feigelson et al. 2005). The scale factors between the 
COUP and CCCP XLFs are 25 and 20, for the lightly 
and heavily obscured samples, respectively, giving a total 
of 45 x 839 = 3.8 x 10 4 diskless CCCP cluster members 
with masses >0.1 M©. This result is a lower limit due to 
the assumption that sources for which Lx could not be 
obtained actually have L x values below the XLF turnoffs 
at log (Lx/[erg s^ 1 ]) ~ 30.6 and 31.6 for the lightly and 
heavily obscured samples, respectively (Figure 11). 

Combining >3.8 x 10 4 diskless stars traced by the XLF 
(Figure 11) with >2 x 10 4 disk-bearing YSOs predicted 
by the YMF (Figure 3) gives >5.8 x 10 4 young stars for 
the global Carina Nebula population, in agreement with 
the estimates of Fll and Smith & Brooks (2007). The 
total stellar mass is 4.9-8.6 x 10 4 M© (the upper bound 
is computed using the XLF scaling analysis from Fll). 
Because the large majority of CCCP sources classified as 
likely Carina members (Broos et al. 2011b) are associated 
with the clustered stellar population (Fll), we deduce 
that the age spread represented by the global XLF is 



similar to the ~5 Myr age spreads measured for the ma- 
jor ionizing clusters (Wang et al. 2011; Wolk et al. 2011, 
and references therein). The time-averaged global SFR 
has been 0.010-0.017 M© yr" 1 over the past -5 Myr. 

The combined YMF and XLF results above lead 
straightforwardly to a global mid-IR excess fraction of 
^30% among the Carina Nebula population. Compared 
to the commonly-used i^g-excess fraction, this mid-IR 
excess fraction is much closer to the astrophysical frac- 
tion of stars with inner disks. YSOs dominated by cir- 
cumstellar disks or envelopes do not separate cleanly 
from reddened, diskless stars in JHKs color-color di- 
agrams (Whitney et al. 2003, RW06). Applying the 
i^g-excess selection criteria of Preibisch et al. (2011) 
to the PCYC sources detected in 2MASS JHK S (Ta- 
ble 2), we find that only 36% of mid-IR selected YSOs 
exhibit a K$ excess. We thus predict a global K$ ex- 
cess fraction among the Carina Nebula population of 
30% x 36% = 11%. This result agrees very well with the 
<10% global ifg-excess fraction reported by Preibisch 
et al. (2011) from deep JHK S imaging of the CCCP 
stellar population, considering that the X-ray detection 
frequency is expected to be lower for YSOs than for disk- 
less PMS stars. We note that the stellar population in 
the Carina Nebula is far from coeval, hence large varia- 
tions from the ~30% global average mid-IR excess frac- 
tion are expected among different sub-populations, for 
example the relatively old Tr 15 cluster (Feinstein et al. 
1980) has a very low disk fraction (see Figure 2 and also 
Wang et al. 2011) while the very young Treasure Chest 
cluster may have a disk fraction near 100% (Smith et al. 
2005). The disk fraction also decreases as a function of 
increasing stellar mass (PW10); the global value is dom- 
inated by the numerous low-mass stars, while the OB 
stellar population generally exhibits a low frequency of 
marginal mid-IR excess emission that likely is not indica- 
tive of circumstellar disks (Pll). 

In the South Pillars region, SP10 observed a general 
pattern of YSOs distributed spatially in the cavities be- 
tween the ends of pillars and the ionizing stars and con- 
cluded that when star formation is triggered by feed- 
back from massive stars, YSOs are left in the wake of re- 
treating, evaporating pillars and molecular cloud surfaces 
eroded by advancing ionization fronts. Just over 50% of 
the YSOs in the PCYC are located in the South Pillars 
field presented by SP10 (Figures 1 and 2). The wide-field 
coverage of the Vela-Carina survey allows us to gener- 
alize this picture of ongoing, potentially feedback-driven 
star formation beyond the South Pillars to the entire 
Carina Nebula. If triggering were indeed the dominant 
mode for ongoing star formation, then the high present- 
day SFR in the Carina Nebula would demonstrate that 
feedback can provide the mechanism for sustaining and 
propagating star formation throughout a CMC complex. 

We caution, however, that while most of the active star 
formation has migrated outward (Smith & Brooks 2007), 
making a clear division between a triggered YSO popula- 
tion and the OB stars that provide the trigger would be 
an oversimplification. The large concentrations of YSOs 
associated with the central, massive Tr 16 and Tr 14 
clusters (Figures 1 and 2) indicate that these clusters 
are hotbeds of recent star formation; they arc still being 
built. In contrast, the older Tr 15 and Bochum 11 clus- 
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ters, both located in the outer regions of the Carina Neb- 
ula, have no associated YSO sub-clusters even though 
both are obvious X-ray star clusters (Figure 2). The 
Carina Nebula is famous for its evolved super-massive 
stellar population, n Car being the capstone, but the 
numerous OB stars apparently associated with compact 
clusters or groups of YSOs (Figure 10, and see also Fig- 
ure 5 of Pll) could be as young as a few 10 5 yr, especially 
considering that the intense winds and radiation fields of 
massive stars likely destroy their natal envelopes over 
very short timescales compared to the disk lifetimes of 
lower-mass stars (Hollcnbach et al. 1994). 

The triggered star formation hypothesis is very difficult 
to test in a region as complicated as Carina; proof would 
require additional information beyond the observed cor- 
relation between YSOs and ionization fronts to estab- 
lish a causal relationship that the ongoing star forma- 
tion could not have occurred in the absence of a trig- 
gering mechanism provided by the massive stars. High- 
resolution intcrfcrometric radio or sub-mm detections of 
dense, massive, molecular cores that lack active star for- 
mation and remain protected from external radiation 
fields within the Carina molecular clouds could provide 
such evidence. 

An alternate hypothesis of turbulence-driven, hierar- 
chical fragmentation and collapse of the Carina CMC 
complex (Sabbi et al. 2007, and references therein) also 
appears consistent with currently available data. This 
idea could explain the cluster-of-clusters morphology of 
the Carina stellar populations and is supported by our re- 
sults that a large, approximately coeval YSO population 
is widely distributed throughout the Carina complex and 
the present-day SFR is approximately the same as the 
time-averaged historical SFR. In this scenario, feedback 
from massive stars need not trigger star formation, but 
it could regulate star formation by rapidly sweeping the 
natal molecular cloud fragments away from young em- 
bedded clusters, revealing and even unbinding these YSO 
clusters (as suggested by SP10) over shorter timescales 
than expected from observations of more pristine molec- 
ular cloud environments such as M17 SWex (PW10). 

In summary, the global star formation history of the 
Carina GMC complex is characterized by an approxi- 
mately constant SFR of 0.010-0.017 M yr^ 1 over the 
past >5 Myr, perhaps set in motion by the hierarchical 
collapse and fragmentation of the GMC into a cluster 
of clusters. This average SFR has been punctuated by 
intense bursts of massive star formation. The largest 
burst produced the super-massive stellar population in 
Tr 16, including n Car, about 3-4 Myr ago (e.g. DeGioia- 
Eastwood et al. 2001; Wolk et al. 2011), followed by a 
second burst that formed the very massive stellar pop- 
ulation in the core of the Tr 14 cluster, perhaps as re- 
cently as 0.5 Myr ago (Sana et al. 2010). Subsequent to 
the formation of the very massive ionizing clusters, the 
recent star formation history in the Carina Nebula has 
been driven or at least regulated by feedback from the 
massive stars. 

7.2. The Production of X-rays by Intermediate-Mass 
PMS Stars 

Our comparative study in §5 between the subpopu- 
lations of the PCYC with and without X-ray matches 
in the CCCP catalog (Broos et al. 2011a) reveals some 



qualitative trends that provide clues into the origin of X- 
ray emission in intermediate- mass PMS stars. Across the 
<1— 10 M Q mass range sampled by the PCYC (Figure 5), 
X-ray emission appears to be activated during the pro- 
tostcllar evolutionary phase, and this emission becomes 
more readily detectable in relatively evolved Stage 0/1 
objects for which the absorption to the central star is 
reduced. It is well-established that the X-ray luminos- 
ity of low-mass, T Tauri stars increases as they evolve 
from disk-dominated (Stage/Class II) to diskless (Class 
III) PMS stars (Stassun et al. 2004; Telleschi et al. 2007; 
Prisinzano et al. 2008). In stark contrast, we find that 
intermediate-mass, disk-bearing YSOs continue to pro- 
duce X-rays as they evolve beyond the Stage 0/1 phase, 
but this emission fades away as the circumstellar disks 
evolve and the central stars pass through a Herbig Ae/Be 
phase during their approach to the main sequence. 

X-ray detections of intermediate-mass, disk-dominated 
YSOs in the mass range 1.6-4 M Q are correlated with 
higher disk mass only for stars with cool photospheric 
temperatures (Figure 7). Although these stars will ar- 
rive on the main sequence as A and B stars with T e s > 
7500 K and fully radiative envelopes, currently they have 
T c ff < 5500 K and spectral types later than G2. Dur- 
ing their early PMS evolution, intermediate-mass stars 
can thus be regarded, in terms of their stellar struc- 
ture and surface temperatures, as more luminous analogs 
of classical T Tauri stars, complete with convective en- 
velopes. We therefore suggest that X-ray emission from 
intermediate-mass YSOs originates in a scaled-up ver- 
sion of the mechanism powering X-rays from T Tauri 
stars; convection-driven magnetic reconnection activity 
(Stassun et al. 2004; Preibisch et al. 2005). Because 
intermediate-mass stars reach the main sequence within 
a few Myr (Siess et al. 2000), this emission mechanism 
is available for only a short time compared to low-mass 
PMS stars. The rapid fading of X-ray emission coincides 
with the rapid dissipation of the circumstellar disk as the 
temperature of the central star rises. 

This rapidly fading dynamo hypothesis can also ex- 
plain a puzzling feature of the CCCP XLF. The slope 
of the lightly-obscured CCCP XLF, fitted with a power- 
law function using the method of Maschbergcr & Kroupa 
(2009), is Txlf = 1-3, significantly steeper than the 
slopes of both the the comparison lightly-obscured 
COUP XLF (r XL F = 0.9) and the heavily-obscured 
CCCP XLF (r XL F = 1-0; Figure 11). If the lightly- 
obscured CCCP XLF traces an older stellar popula- 
tion compared to either the COUP or the heavily- 
obscured CCCP XLFs, then a relatively rapid decay in 
Lx among the more luminous, intermediate-mass stars 
would steepen the observed XLF slope, just as mass- 
dependent circumstellar disk evolution steepens the ob- 
served YMF slope (Figure 3). A similar deficit of stars 
with \ogLx > 31 erg s _1 steepens the observed XLF 
slope of the Tr 15 cluster (Wang et al. 2011). With an 
age >5 Myr, Tr 15 represents a significantly older stel- 
lar population compared to cither the ONC or the other 
massive clusters in the Carina Nebula. 

8. CONCLUSIONS 

We have used the Spitzer Vela-Carina point-source 
lists, in combination with complementary IR datasets, 
to carry out an unbiased search for YSOs with IR excess 
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emission throughout the Carina Nebula. We produce a 
Pan-Carina YSO Catalog (PCYC) of 1439 highly-reliable 
candidate YSOs, 410 of which have X-ray counterparts 
from the CCCP catalog. The PCYC comprises one of the 
largest samples of intermediate-mass YSOs in a single 
massive star-forming complex compiled to date, and it is 
the largest such collection with a complete set of comple- 
mentary X-ray observations. From our initial analysis of 
this rich dataset, we can draw the following conclusions: 

• The mass function (YMF) of the Carina YSOs ex- 
hibits a power-law slope Tymf = 3.3, significantly 
steeper than the Salpeter-Kroupa IMF. Such a 
steep slope, first observed for the YSO population 
of the M17 SWex IR dark cloud, is at least partially 
the result of an observational selection effect caused 
by rapid disk evolution among intermediate-mass 
YSOs (PW10). M17 SWex and Carina are very 
different environments; the former is a cold, IR 
dark cloud apparently experiencing its first wave 
of star formation while the latter is an evolved 
H II region powered by some of the most massive 
stars in the Galaxy. While the YMF slopes are sim- 
ilar, the lower critical mass of the PCYC sample, 
below which the YMF may agree with a normal 
Salpeter-Kroupa IMF, likely reflects a larger age 
spread among the YSOs in the Carina Nebula. 

• The time-averaged global SFR in the Carina com- 
plex has been 0.010-0.017 Msun yr" 1 over the past 
<5 Myr, punctuated by intense bursts of activity 
that produced the very massive stellar populations 
in the Tr 16 and Tr 14 clusters. If the present- 
day Galactic SFR is a few M yr _1 , then the star 
formation activity of the entire Milky Way is equiv- 
alent to a few hundred Carina complexes. The 
cluster-of-clusters pattern of star formation in the 
Carina Nebula is consistent with hierarchical frag- 
mentation and collapse of a GMC complex (Sabbi 
et al. 2007). 

• The present-day global star formation rate in the 
Carina Complex is >0.008 M Q yr _1 . The mas- 
sive central ionizing clusters Tr 16 and Tr 14 still 
host intense concentrations of ongoing star for- 
mation. Ongoing star formation activity is dis- 
tributed widely throughout Carina, driven or regu- 
lated by feedback from the massive stellar popula- 
tion, which acts to quickly reveal and even unbind 
young, embedded clusters (SP10). 

• The PCYC includes 62 candidate X-ray-emitting 
Stage 0/1 YSOs, a sample large enough to provide 
statistical evidence that X-ray detections correlate 
with more evolved protostars. Our results support 
previous findings that the physical mechanism pro- 
ducing X-rays in PMS stars is activated during the 
protostellar evolutionary phase, and circumstellar 
material contributes significantly to the absorption 
of X-rays from the youngest YSOs (Getman et al. 
2007; Prisinzano et al. 2008). 

• The apparent X-ray emission from intermediate- 
mass YSOs is usually consistent with the presence 
of unresolved, low-mass companions. But we also 



isolate a class of intermediate-mass Stage II YSOs 
with cooler photospheres and relatively massive cir- 
cumstellar disks for which the X-ray emission ap- 
pears to be intrinsic. We suggest that the progen- 
itors of Hcrbig Ae/Be stars produce X-rays dur- 
ing their early PMS evolution, perhaps powered by 
magnetic reconnection activity during a convective 
atmosphere phase, but this mechanism decays as 
the stars approach the main sequence. A rapid fad- 
ing of the stellar dynamo could thus act to steepen 
observed XLF slopes as a function of increasing 
cluster age (Wang et al. 2011). 

• We identify 3 new very young, tight YSO sub- 
clusters in the northwestern region of the Carina 
Nebula, 2 of which are associated with significant 
concentrations of X-ray sources (Fll). 

The full PCYC dataset is available in the electronic 
versions of Tables 2 and 3. We hope that this catalog 
will serve as a valuable database for future multiwave- 
length observations and modeling of star formation in 
the Carina Complex. 
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Pan-Carina YSO Catalog: Basic IR Observational Data a 
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Note. — This Tabic is available in its entirety in a machine- read able form in the online journal. A portion is reproduced here for guidance 
regarding its form and content and to highlight sources of interest. Sources are grouped according to their apparent association with the YSO 
sub-clusters discussed in §6. 

a The typical photometric uncertainty is 0.1 mag. 

k Source names arc preceded by cither a "C" for the highly-reliable Vela-Carina Catalog or an "A" for the morc-complctc Archive. 

C Quality flags are defined as follows: — non-detection, 1 — detection, 2 — flux density lower limit (magnitude upper limit), and 3 — flux density 

upper limit (magnitude lower limit). 
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(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


PCYC 


\ lu 6 -^bol/ 




(M+) 


<j(M+) 


(lo 


(J M.„,,\ 


cr(\c>& fVf^T,^r \ 

V 6 cnv } 




CCCP Match 


Membership* 




1NO. 


/T ~\ 
( L ©) 


( L o) 


/TV T \ 

(M ) 


( M ©J 




3 yi J 


(M yr ) 


btage 




Cl/Gr 


Reg 


INotcs 


642 


1.9 


0.8 


2.8 


0.9 








A 




C5 


A 


Cr 232 


654 


1.8 


1.1 


2.2 


1.5 




-5.2 


2.0 


0/1 






A 


Cr 232 


655 


2.1 


0.3 


3.9 


1.3 




-5.0 


2.1 


A 






A 


Cr 232 


656 


2.3 


1.4 


3.1 


1.8 




-5.2 


1.9 


A 






A 


Cr 232 


657 


2.2 


1.3 


3.3 


1.6 




-5.3 


1.9 


A 






A 


Cr 232 


660 


1.8 


1.1 


2.2 


1.2 




-8.5 


2.4 


A 




C5 


A 


Cr 232 


699 


1.0 


0.8 


0.7 


1.0 




-4.8 


2.2 


0/1 


104451.91-602511.9 


G18 


D 


3X 


1295 


1.7 


0.1 


3.4 


0.2 








III 






D 


o 


1297 


1.2 


0.4 


1.3 


1.3 




-5.0 


2.3 


0/1 






D 


o 


1302 


1.4 


1.1 


1.9 


0.6 








A 






D 


O 


1304 


2.6 


1.5 


4.2 


1.5 








A 






D 


O 


1308 


2.0 


1.2 


3.1 


0.9 








II 






D 


O 


1312 


2.2 


0.1 


3.7 


0.4 








11 


104748.69-602611.1 




D 


O 


1313 


1.3 


0.8 


2.8 


0.6 








II 






D 


o 


1315 


1.5 


0.4 


2.1 


1.2 




-5.0 


1.9 


0/1 


104750.03-602619.1 




D 


o 


1316 


1.2 


0.8 


1.8 


1.1 




-6.3 


2.5 


A 






D 


o 


1317 


2.0 


1.2 


1.4 


1.5 




-4.6 


1.9 


0/1 






D 


o 


1319 


1.2 


0.5 


1.4 


1.0 




-4.9 


1.9 


0/1 






D 


o 


1320 


2.6 


1.5 


4.3 


1.5 




-3.3 


1.5 


0/1 






D 


o 



Note. — This Tabic is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding 
its form and content and to highlight the sources of interest presented in Table 2 and §6. 

a The Membership columns follow Fll. 

k The Notes column includes the following: O, P, Q, R, or Cr 232 refer to membership in the groups discussed in §6; 2X or 3X indicates that the 
Spitzer source was matched to 2 or 3 CCCP X-ray sources; and MSX means that the YSO is also an MSX source. 

Table 4 

Candidate X-ray- Emitting Stage 0/1 YSOs 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 

Physical Parameters from YSO Models Observed X-ray Properties' 3 



PCYC <M*) a(M t ) (logT*> (r(logT*) (log ^) (A v , t ) Max,<r(A Vi i) a CCCP Match Net B med ian logF t 
No. (M ) (M Q ) (K) (K) (yr- 1 *) (mag) (mag) CXOGNC J. Counts (keV) (erg/cm 2 /s) 



Hard X-ray Counterparts (-E mo dian > 2 keV) 



138 


3.1 


0.7 


3.6 


0.0 


-5.0 


5.0 x 10 2 


4.5 x 10 2 


104248.10-592539.3 


101.5 


3.3 


-13.4 


140 


0.8 


0.9 


3.5 


0.1 


-5.8 


4.1 x 10° 


6.4 x 10 2 


104248.50-594835.5 


5.2 


2.3 


-14.9 


179 


10.1 


0.9 


3.6 


0.0 


-3.6 


7.2 x 10 2 


4.9 x 10 2 


104309.83-592454.2 


44.9 


4.0 


-13.9 


239 


1.2 


1.2 


3.6 


0.2 


-5.6 


4.6 x 10 1 


4.7 x 10 2 


104326.70-602622.7 


149.6 


4.7 


-12.9 


282 


1.9 


1.5 


3.7 


0.2 


-5.6 


3.1 x 10 1 


9.9 x 10 4 


104339.83-593300.9 


3.6 


3.0 


-14.9 


318 


2.8 


1.0 


3.7 


0.2 


-5.8 


1.5 x 10 1 


1.5 x 10 1 


104347.79-593304.4 


5.6 


2.9 


-14.8 


348 


3.7 


2.4 


3.8 


0.4 


-5.4 


4.9 x 10 1 


4.0 x 10 5 


104351.67-593257.0 


9.6 


2.1 


-14.4 


353 


1.9 


1.1 


3.7 


0.1 


-5.8 


7.0 x 10° 


5.1 x 10 3 


104352.89-593921.8 


92.9 


2.4 


-13.6 


371 


4.0 


1.1 


3.7 


0.1 


-5.4 


2.5 x 10 1 


1.7 x 10 1 


104355.73-592417.5 


55.3 


2.8 


-13.6 


373 


3.7 


1.2 


3.7 


0.0 


-5.3 


2.9 x 10 1 


1.3 x 10 1 


104355.80-592421.7 


11.8 


3.2 


-14.1 


387 


1.5 


1.7 


3.6 


0.2 


-5.4 


5.7 x 10 1 


1.7 x 10 3 


104357.77-592833.3 


8.5 


2.2 


-14.6 


419 


2.8 


1.6 


3.6 


0.1 


-4.8 


1.8 x 10 2 


1.6 x 10 s 


104403.08-593748.7 


19.0 


2.6 


-14.5 


484 


9.5 


0.8 


4.0 


0.6 


-5.2 


4.1 x 10 1 


3.5 x 10 1 


104417.86-602745.9 


5.5 


4.7 


-14.4 


598 


2.2 


1.6 


3.7 


0.3 


-5.3 


3.7 x 10 1 


1.0 x 10 3 


104437.14-595802.4 


29.1 


4.0 


-13.9 


623 


3.5 


1.3 


3.8 


0.5 


-4.9 


4.1 x 10 1 


9.3 x 10 1 


104440.31-594339.7 


2.3 


5.4 


-14.9 


686 


0.9 


1.0 


3.6 


0.1 




5.2 x 10° 


1.5 x 10 2 


104449.16-595346.6 


138.3 


2.9 


-13.5 


699 


0.7 


1.0 


3.5 


0.1 


-4.7 


3.0 x 10 1 


1.4 x 10 3 


104451.91-602511.9 


376.8 


5.2 


-12.3 


713 


1.6 


0.7 


3.7 


0.1 


-6.2 


8.9 x 10° 


7.0 x 10° 


104455.40-592634.5 


35.5 


4.4 


-14.1 


774 


2.1 


1.1 


3.7 


0.1 


-6.6 


1.5 x 10 1 


1.6 x 10 1 


104507.15-600156.9 


79.2 


3.0 


-13.8 


897 


0.6 


0.5 


3.6 


0.1 


-5.3 


3.2 x 10 1 


7.3 x 10 2 


104530.83-595809.5 


5.3 


2.2 


-15.0 


1005 


2.4 


0.9 


3.7 


0.1 


-9.0 


1.5 x 10 1 


8.9 x 10° 


104549.19-601637.3 


19.4 


3.0 


-14.3 


1073 


1.7 


1.0 


3.6 


0.0 


-5.7 


8.4 x 10° 


4.9 x 10° 


104603.53-595339.3 


28.6 


2.2 


-14.3 


1185 


1.0 


0.9 


3.7 


0.2 


-5.8 


1.5 x 10 1 


2.6 x 10 1 


104637.74-594804.4 


15.9 


2.2 


-14.6 


1236 


2.7 


0.4 


3.6 


0.0 


-5.0 


2.8 x 10 1 


6.1 x 10° 


104659.85-602535.0 


10.3 


2.3 


-14.4 


1315 


2.1 


1.2 


3.6 


0.0 


-5.3 


2.7 x 10 1 


1.6 x 10 1 


104750.03-602619.1 


17.2 


2.7 


-14.0 


1367 


7.6 


2.5 


4.0 


0.2 


-5.4 


5.6 x 10 2 


1.7 x 10 6 


104824.37-600759.5 


22.0 


4.2 


-13.9 


Soft X-ray Counterparts (i? mo( jian < 2 keV) 


83 


0.8 


0.6 


3.6 


0.0 


-5.1 


2.1 x 10 1 


2.9 x 10 1 


104216.00-593607.6 


13.5 


1.5 


-14.4 


148 


0.6 


0.5 


3.6 


0.1 


-4.5 


4.6 x 10 1 


1.4 x 10 2 


104251.41-592535.6 


5.5 


1.5 


-14.9 
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(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 

Physical Parameters from YSO Models Observed X-ray Properties' 3 



PCYC (M*> CT (M*) (logT*) a(logT*) (log ^> (A v , t > Max,<r(Ay, t ) a CCCP Match Net £ mcdi an logF t 



No. 


(M ) 


(M ) 


(K) 


(K) 


(yr- 1 ) 


(mag) 


(mag) 


CXOGNC J. Counts 


(keV) 


(erg/cm- 


lOU 


1 n 

1 .U 


n 7 

u. / 


O. 1 


n 9 

U.Z 


^ A 


2.8 X 


10 1 


2.9 x 


10 2 


104255.40-592456.5 


3.6 


i .o 


-14 9 




1 2 


n 7 

u. / 


O.U 


U. 1 


-u.o 


1.0 x 


10 1 


3.3 x 


10 1 


104312.59-592516.2 


30.1 


1 .o 


-14 4 


9&& 
Zoo 


9 

z.o 


i .u 


o.u 


n n 
u.u 


-o. i 


2.5 x 


10 2 


5.4 x 


10 3 


104341.07-594537.2 


51.1 


i .o 


-14 4 


346 


2.5 


0.6 


3.6 


0.0 


-4.7 


3.1 x 


10 2 


6.6 x 


10 2 


104351.59-600525.2 


17.0 


1.3 


-14.6 


352 


1.7 


1.6 


3.6 


0.2 


-5.7 


2.3 x 


10 1 


6.5 x 


10 4 


104352.58-593025.1 


8.7 


1.4 


-15.1 


OOO 


9 n 

z.u 


n q 
u.y 


3 7 

O. 1 


U. 1 


-o.o 


3.2 x 


10 1 


5.0 x 


10 4 


104353.56-593410.1 


5.6 


i .o 


-lO.U 




2 1 


1 ^ 
1 .0 


3 7 

O. 1 


U. 1 


-u.o 


4.4 X 


10° 


1.0 x 


10 3 


104355.20-600519.8 


11.5 


1 2 


-14 4 


o / u 


u.o 


7 
u. / 


o.u 


u.o 




4.8 x 


10° 


4.2 x 


10 1 


104355.92-592537.6 


19.5 


i .o 


14 7 
-14. ( 


424 


2.8 


1.1 


3.7 


0.2 


-5.5 


5.6 x 


10° 


5.5 x 


10° 


104403.87-593308.3 


4.7 


1.0 


-15.4 


432 


2.3 


0.9 


3.7 


0.0 


-5.3 


1.2 x 


10 2 


2.0 x 


10 2 


104406.04-595510.6 


11.1 


1.7 


-14.6 




9 3 

z.o 


1 2 


3 7 

O. 1 


n 9 
u.z 


^ 7 


2.8 x 


10° 


1.8 x 


10 1 


104409.95-593436.1 


27.2 


l.O 


AAL 

- 14. 


^oy 


1 .O 


n 7 

u. / 


o.u 


U. 1 


7 9 
- / .z 


8.7 x 


10° 


5.4 x 


10 1 


104411.75-592431.1 


19.5 


1 .u 


-14. O 


4yu 


1 7 
1 . 1 


1 4 


o.u 


U. 1 


-o.o 


4.5 x 


10 1 


1.6 x 


10 2 


104421.14-593000.9 


4.6 


l.O 


-±o. 1 




U.O 


1 o 
i .u 


o.u 


U. 1 


9 

-u.z 


3.4 x 


10° 


3.1 x 


10 2 


104423.08-600819.3 


26.9 


1 2 


-14 4 


508 


0.7 


0.9 


3.6 


0.2 


-5.8 


5.8 x 


10° 


6.3 x 


10 2 


104423.32-592513.5 


5.5 


1.6 


-15.1 


521 


1.3 


1.1 


3.7 


0.3 


-8.3 


5.1 x 


10° 


2.0 x 


10 3 


104425.68-593923.5 


55.1 


1.3 


-14.5 


040 


1 7 
1 . 1 


1 7 


o.u 


U. 1 


-0.4 


5.7 x 


10 1 


6.7 x 


10 2 


104429.88-594533.8 


4.6 


l.O 


-15 2 




q n 
y .u 


^ 9 
O.z 


4.u 


u.u 


A Q 


8.7 x 


10 1 


2.4 x 


10 3 


104431.23-593311.2 


18.2 


1 Q 

i .y 


1 A ^ 

-14. 




1 ^ 

l.O 


1 9 
i .z 


3 7 

O. 1 


n 9 

u.z 


7 7 


5.6 x 


10° 


1.1 X 


10 4 


104436.06-594258.4 


12.9 


l.O 


1 A 7 
-14. ( 




1 ^ 
1 .O 


u.o 


o.u 


n n 
u.u 


^ A 

-0.4 


3.7 x 


10 1 


1.3 x 


10 3 


104446.07-595816.8 


12.6 


l.O 


1/1 Q 


694 


1.5 


1.5 


3.6 


0.2 


-5.5 


4.0 x 


10 1 


1.3 x 


10 s 


104450.81-594442.9 


58.2 


1.6 


-14.2 


768 


2.6 


1.4 


3.7 


0.1 


-5.9 


1.1 X 


10 2 


4.6 x 


10 2 


104505.75-594405.0 


16.4 


1.3 


-14.7 


795 


1.3 


0.9 


3.7 


0.2 


-6.2 


4.3 x 


10° 


2.4 X 


10° 


104511.09-594533.5 


66.4 


1.4 


-14.1 


890 


0.7 


0.9 


3.6 


0.2 


-5.8 


5.9 x 


10° 


1.4 x 


10 3 


104529.17-601357.0 


20.7 


1.4 


-14.5 


928 


1.0 


0.8 


3.6 


0.1 


-6.0 


9.7 x 


10° 


3.0 x 


10 2 


104536.94-595611.3 


9.5 


2.0 


-14.6 


987 


0.6 


0.7 


3.5 


0.2 


-5.3 


3.3 x 


10 1 


9.1 x 


10 2 


104546.73-600343.3 


14.7 


0.9 


-15.1 


1025 c 


4.2 


0.5 


3.7 


0.0 


-5.3 


5.1 x 


10 1 


6.2 x 


10 1 


104553.71-595703.9 


78.1 


1.5 


-13.8 


1055 


3.6 


2.8 


4.0 


0.6 


-5.2 


7.7 x 


10° 


5.1 x 


10° 


104559.44-600517.9 352.0 


1.6 


-13.2 


1141 


2.5 


0.8 


3.7 


0.1 


-5.8 


1.2 x 


10 1 


1.8 x 


10 1 


104622.54-594523.9 


14.2 


1.7 


-14.3 


1193 


0.5 


0.7 


3.5 


0.2 


-5.8 


4.0 x 


10° 


1.0 x 


10 3 


104642.52-600826.7 


7.9 


1.5 


-14.9 


1204 


1.3 


0.7 


3.6 


0.1 


-5.5 


4.4 X 


10 1 


1.8 x 


10 3 


104646.68-600932.3 


4.0 


1.4 


-14.9 


1276 


0.8 


0.8 


3.6 


0.1 


-8.1 


7.7 x 


10° 


2.8 x 


10 1 


104717.42-600906.1 


61.0 


1.6 


-13.9 


1390 


0.4 


0.6 


3.5 


0.1 


-5.9 


4.0 x 


10° 


2.6 x 


10 2 


104836.21-601408.7 


7.7 


1.3 


-14.9 



a The probability distributions of Av,t can be strongly skewed, with long tails to high values. This often yields very large values for a(Ay,t), 
hence in cases where <r(Av,t) ^> (Av.t) we define M.ax(Av,t) = &{Av,t)- 

k Columns (10)— (12) list the quantities NetCounts t , EncrgyFlux t , and MedianEnergy t from Table 1 of Broos ct al. (2011a). 

C PCYC 1025 is located in the the very crowded center of the Treasure Chest, hence it is very likely a combination of multiple confused sources 
in the IRAC beam, and the physical association of the CCCP source with the IRAC source is highly questionable. 



